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Abstract 

a.: 

Recent achievements in the heavy quark theory are critically reviewed. The emphasis 
is put on those aspects which either did not attract enough attention or cause 
^ ■ heated debates in the current literature. Among other topics we discuss (i) basic 

parameters of the heavy quark theory; (ii) a class of exact QCD inequalities; (iii) 
new heavy quark sum rules; (iv) virial theorem; (v) applications (|V c t,| from the 
total semileptonic width and from the B — > D* transition at zero recoil). In some 
instances new derivations of the previously known results are given, or new aspects 
addressed. In particular, we dwell on the exact QCD inequalities. Furthermore, a 
toy model is considered that may shed light on the controversy regarding the value 
of the kinetic energy of heavy quarks obtained by different methods. 



* Permanent address: Petersburg Nuclear Physics Institute, Gatchina, 
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1 Introduction 



Quark-gluon dynamics are governed by the QCD Lagrangian 

C = ~\g%G% + J> #g + £ Q(i V> - m Q )Q = 

1 Q 

Aight + E QVP - ™q)Q 

Q 

where is the gluon field strength tensor, the light quark fields (u, d and s) are 
generically denoted by q and are assumed, for simplicity, to be massless, and the 
heavy quark fields are generically denoted by Q. To qualify as a heavy quark Q 
the corresponding mass term tuq must be much larger than Aqcd- The charmed 
quark c can be called heavy only with some reservations and, in discussing the heavy 
quark theory, it is more appropriate to keep in mind h quarks. The hadrons to be 
considered are composed from one heavy quark Q, a light antiquark q, or diquark 
qq, and a gluon cloud which also contains light quark-antiquark pairs. The role of 
the cloud is, of course, to keep all these objects together, in a colorless bound state 
which will be generically denoted by Hq. 

The light component of Hq, its light cloudfj has a complicated structure - the 
soft modes of the light fields are strongly coupled and strongly fluctuate. Basi- 
cally, the only fact which we know for sure is that the light cloud is indeed light; 
typical frequencies are of order of Aq CD . One can try to visualize the light cloud 
as a soft medium. The heavy quark Q is then submerged in this medium. If the 
hard gluon exchanges are discarded, the momentum which the heavy quark can 
borrow from the light cloud is of order of Aqcd- Since it is much smaller than 
rriQ, the heavy quark-antiquark pairs can not play a role. In other words, the field- 
theoretic (second-quantized) description of the heavy quark becomes redundant, and 
under the circumstances it is perfectly sufficient to treat the heavy quark quantum- 
mechanically. This is clearly infinitely simpler than any field theory. Moreover, 
one can systematically expand in I/itlq. Thus, in the limit rriQ /Aqcd ~~ * oo the 
heavy quark component of Hq becomes easily manageable, allowing one to use the 
heavy quark as a probe of the light cloud dynamics. Treating the heavy quark Q 
in Hq as a non-relativistic object submerged in a soft gluon background field of the 
light cloud, we open the door for the use of a large variety of methods developed 
in Quantum Electrodynamics long ago [[l], |J, for instance, the Pauli expansion, the 
Foldy-Wouthuysen technique and so on. 

The special advantages of the limit rriQ — > oo in QCD were first emphasized by 
Shuryak . The next logical step was the observation of the heavy quark symmetry 
0, |5|, |J. The heavy quark theory in QCD was finally formalized in Refs. where 

^dn some papers devoted to the subject the light cloud is referred to as "brown muck" , which 
seems to be an unfair name for the soft components of the quark and gluon fields - perhaps we 
have not yet unveiled their beauty. 
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the systematic I/tuq expansion of the Hq dynamics was cast in the form of the 
effective Largangians. Applications of the heavy quark theory are so numerous 
that they comprise now a significant and, perhaps, the most active part of QCD, 
with a well-known record of recent successes and breakthroughs. Many exhaustive 
reviews are devoted to the subject ||, and we feel there is no need to repeat what 
has already become common knowledge in our community. Our task, as we see 
it, is to concentrate on several aspects of the heavy quark theory which are still 
controversial. These aspects, as a rule, are related to a specific structure of QCD 
where the gluon field has two faces. The gluon degrees of freedom play a role of a soft 
background medium, and, simultaneously, they are responsible for the radiative a s 
corrections reflecting the presence of the hard component in the gluon field. These 
two components - soft and hard - are tightly intertwined, and in many instances it 
is hard to disentangle them. Without untangling the hard component one can not 
build consistent 1 / rriQ and a s expansions. This is in sharp distinction with, say, the 
Pauli expansion in QED where this problem simply does not arise. The background 
electromagnetic field can be taken as soft as we want, the radiative corrections are 
fully calculable and small (i.e. they can be neglected in most of the problems). 

The most remarkable example of the problem mentioned above is the definition 
of the heavy quark mass - what quantity appears in the 1/ttiq expansion after all? 
In electrodynamics the electron mass is a quantity that can be measured. Even in 
this case, though, one encounters some (mild) problems in theoretical description, 
due to the masslessness of photons. As a matter of fact, the electron never appears 
"alone" ; it is always surrounded by a photonic cloud, manifesting itself in infrared 
divergences. The strategy allowing one to treat these divergences and define the 
mass of an "isolated" electron is well-known. The situation becomes much more 
complicated in QCD, where isolated quarks simply do not exist. In defining the 
heavy quark mass one can peel off the gluonic cloud only up to a point. Namely, 
we can eliminate from the definition the contribution coming from hard gluons. 
The contribution associated with soft gluons has to be included in the definition 
of the heavy quark mass simply because we do not know the precise structure of 
the soft-gluon component. How exactly could one separate the hard and soft-gluon 
components in defining the heavy quark mass? This question causes debate in the 
current literature, and will be one of the focal points of the present review. 

Another key parameter of the heavy quark theory is the kinetic operator. For- 
mally its expectation value can be defined as 

& = 7^—{H Q \Q7f 2 Q\H Q ) , jf=-iD, (1.2) 

ZM H Q 

where D stands for the covariant derivative (note: here and elsewhere, we always 
assume, unless explicitly stated to the contrary, that we are in the rest frame of the 
heavy hadron). 

The physical meaning of /x^ is quite evident: the heavy quark inside Hq expe- 
riences a zitterbewegung due to its coupling to other degrees of freedom inside the 
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hadron. Its average spatial momentum squared is /i^. The kinetic energy opera- 
tor appears in the effective Lagrangian in the next-to-leading order, and basically 
determines the magnitude of I/tuq corrections in a large number of expressions 
of practical interest - from the simplest formula for the heavy hadron masses to 
B — ► £)W formfactor at zero recoil. When one tries to proceed from formal expan- 
sions to practical analysis including perturbative 0(a k s ) effects, the same question 
- what contribution is to be included in jj^ - immediately resurfaces. One can find 
in the literature a spectrum of answers. The numerical value of remains rather 
uncertain. Various determinations of this crucial parameter contradict each other. 
We will present a balanced discussion of competing definitions trying to outline a 
procedure which seems fully consistent in QCD. 

A related topic we would like to consider is a QCD analog of the virial theorem. 
The issue was raised recently by Neubert 0, who noted that ^ is related to a 
certain transition amplitude induced by the chromoelectric field. This observation 
was exploited for an alternative determination of /z^, which brought a surprise. A 
numerical value of /i^ emerging in this way was significantly lower than the results 
of some other determinations. The discrepancy is a subject of ongoing discussions. 
We present a transparent derivation of the virial theorem in QCD, and comment on 
the possible causes of the discrepancy. 

This review is organized as follows. First, a brief survey is given of the operator 
product expansion (OPE) as it is used in the heavy quark theory. This survey 
allows us to introduce all basic tools and formulate theoretical problems to be dealt 
with in the remainder of the review as they appear in their natural environment. 
We then proceed to relatively recent applications of the heavy quark theory in 
the inclusive decays of heavy hadrons. The following topics are discussed in some 
detail: (i) basic parameters of the heavy quark expansion; (ii) a class of exact 
QCD inequalities; (iii) the virial theorem and (iv) applications (\V c b\ from the total 
semileptonic width and from the exclusive B — > D* transition at zero recoil). Rather 
than attempting a comprehensive coverage of a variety of results obtained in the last 
two or three years, we dwell on these selected topics which either have not attracted 
sufficient attention so far, or are subject to ongoing debate. A significant part of 
our presentation is based on the results of Ref. [TD[ that provides a consistent field- 



theoretic framework for treating various aspects of the heavy quark theory that are 
under intense investigation at present. 

Even limiting ourselves to this rather narrow task we had to leave aside, with 
heavy heart, many relevant issues, e.g. duality, inclusive non-leptonic decays of 
heavy flavors, etc. We apologize to the authors whose results, due to space limita- 
tions, are left outside the scope of the present review. The interested readers are 
referred to recent talks [ITT . 



3 



2 Operator Product Expansion; Effective Hamil- 
tonian 



The basic theoretical tool of the heavy quark theory is the Wilson operator product 
expansion [ l2| . It is used in two different directions. First, one can restrict oneself 
to the case of a fixed number of heavy quarks Q since QQ fluctuations can be 
ignored. Most often one considers the sector with a single heavy quark Q which 
is treated as being submerged into a gluonic medium acting as the sole source of 
its interactions. At this stage we disregard any other possible interaction of the 
heavy quark, e.g. electromagnetic, weak and so on. The original QCD Lagrangian 
( p. . 1|) is formulated at very short distances, or, which is the same, at a very high 
normalization point fi = Mo, where Mo is the mass of an ultraviolet regulator; i.e., 
the normalization point is assumed to be much higher than all mass scales in the 
theory, in particular, fi ^> vtlq. An effective theory for describing the low-energy 
properties of heavy flavor hadrons is obtained by evolving the Lagrangian from the 
high scale M down to a normalization point /i lying below the heavy quark masses 
rriQ. This means that we integrate out, step by step, all high-frequency modes in the 
theory thus calculating the Lagrangian The latter is a full-fledged Lagrangian 

with respect to the soft modes with characteristic frequencies less than fi. The 
hard (high-frequency) modes determine the coefficient functions in £(fi), while the 
contribution of the soft modes is hidden in the matrix elements of (an infinite set) 
of operators appearing in C(fi). The value of this approach, outlined by Wilson long 
ago |12| , has become widely recognized and exploited in countless applications. The 
peculiarity of the heavy quark theory lies in the fact that the in and out states contain 
heavy quarks. Although we integrate out the field fluctuations with the frequencies 
down to fi <C rriQ, the heavy quark fields themselves are not integrated out since 
we consider the sector with heavy-flavor charge 7^ 0. The effective Lagrangian C(fi) 
acts in this sector. Since the heavy quarks are neither produced nor annihilated, 
any sector with the given Q-number can be treated separately from all others, as a 



'vacuum" . 



If QCD were solved we could include all modes down to /x = in our explicit 
evaluation of the effective Lagrangian £(//). The resulting Lagrangian would be built 
in terms of the fields of physical mesons and baryons rather than quarks and gluons 
- the latter become irrelevant degrees of freedom in the infrared limit jj, — > 0. All 
conceivable amplitudes could be read off directly from such an effective Lagrangian 
and could be compared with experimental data. 

This picture is of course Utopian: real QCD is not solved in closed form, and 
in explicit calculations of the coefficients in the effective Lagrangian one cannot put 
/i = 0. For decreasing values of /i a larger and larger part of the dynamics has to be 
accounted for in the explicit calculation. We would like to have /1 as low as possible, 
definitely, /i <C niQ. The heavy quark can be treated as a non-relativistic object 
moving in a soft background field only if the latter condition is met. On the other 
hand, to keep computational control over the explicit calculations of the coefficient 



4 



functions we must stop at some \i ^> Aq C d, so that a s (//)/7r is still a sufficiently 
small expansion parameter. In practice this means that the best choice (which we 
will always adopt) is \i ~ several units times Aq CD , i.e. 0.7 to I GeV. All coefficients 
in the effective Lagrangian obtained in this way will be functions of \l. 

Since /i is an auxiliary parameter, predictions for physical quantities cannot 
depend on /i, of course. The \i dependence of the coefficients must be canceled by 
that coming from the physical matrix elements of the operators in C(fi). However, 
in calculating in the hard and soft domains (i.e. above /i and below \i) we make 
different approximations, so that the exact /j independence of the physical quantities 
does not take place. Since the transition from hard to soft physics is very steep, one 
may hope that our predictions will be insensitive to the precise choice of /j provided 
that n ~ several units times Aq CD . 

In descending from Mq to fj,, the form of the Lagrangian ( |1 . 1|) changes, and a 
series of operators of higher dimension appears. For instance, the heavy quark part 
of the Lagrangian takes the form 



£ h eavy = £ { Q{i P ~ ™q)Q + ^Q^G^Q + £ ^QTQqTq \ + O (\ 

1 • 2m Q 2 T\1j m n I V rrir 



(2-1) 

(r) 

where cq and dg^ are coefficient functions, = gG't v t a and t a is the color generator 
(the coupling constant is included into G). We often use the short-hand notation 
oG = o^G^y = ^pTIvG^y. The sum over the light quark flavors q is shown explicitly 
as well as the sum over possible structures T of the four-fermion operators. All 
masses and couplings, as well as the coefficient functions cq and S r \ depend on the 
normalization point. 



The operators of dimension five and higher in Eq. (2.1) are due to the contri- 
bution of hard gluons, with virtual momenta from \i up to Mq. Here the 1/tuq 
expansion is explicit. Implicitly, a 1/rriQ expansion is generated also by the first 
(tree-level) term in the Lagrangian fl2.I|) , 



Cl avy = Q(f-m Q )Q. (2.2) 

Although the field Q in this Lagrangian is normalized at a low point /i, it carries a 
hidden large parameter, mq. Indeed, the interaction of the heavy quark with the 
light degrees of freedom enters through = iD^. The background gluon field 
is weak if measured at the scale mq, which means, of course, that there is a large 
"mechanical" part in the x dependence of Q(x), known from the very beginning, 

Q(x) = e - im ^Q(x) ; (2.3) 

Q(x) is a rescaled bispinor field which, in the leading approximation, carries no 
information about the heavy quark mass. It describes a residual motion of the heavy 
quark inside the heavy hadron [0] with typical momenta of order Aqcd- Remnants 
of the heavy quark mass appear in Q only at the level of I/ttlq corrections. 
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Equation ( |2.3| ) is written in the rest frame of Hq. In an arbitrary frame one 
singles out the factor exp(— irnqv^x^) 



Q{x) 



Q{x) 



(2.4) 



where v M is the four-velocity of the heavy hadron. 

The covariant momentum operator P M acting on the original field Q, when ap- 
plied to the rescaled field Q, is replaced by the operator rriQV^ + ir^, 



iD^Q(x) = e 



mQVf, + iD^j Q(x) = e 



^ (mQV^ + vr^) Q(x) . (2.5) 



If not stated otherwise, we use the rescaled field Q, omitting the tilde in all expres- 
sions where there is no risk of confusion. The rescaled field Q is a four- component 
Dirac bispinor, not a two component non-relativistic spinor which is usually in- 
troduced in the heavy quark effective theory (HQET) 0. The Dirac equation 
(If — 17Iq)Q = can be rewritten as follows in terms of Q: 



7o 



Q 



i 



2m. 



Q 



Q 



7T 2 + \oG 

2mQ 



Q 



(2.6) 



(2.7) 



Armed with this knowledge one can easily expand £° eavy , a ^ the tree level, 
through order 1/m^, 



heavy 



Q 



1 + 7o 



57717 



Q(iip-m Q )Q 

z 

(-(DE) + a- {E X 7? - 7T X E} 



1 + 



Sm 2 Q 



0"7T 



57777 



7T - 



2777Q 



7T(J - 



2 \ 777?, . 



(2- 



where a denote the Pauli matrices and (vr<?) 2 = if 2 + oB , E and .B denote the 
background chromoelectric and chromomagnetic fields, respectively, with the cou- 
pling constant g and the color matrix t a included in the definition of these fields. 
There is nothing new in this Lagrangian; at the tree level it is the same as in QED. 
The non-relativistic expansions in QED have been known since the thirties, see e.g. 
Chapter 4 of Bjorken and Drell U or Sect. 33 of the Landau and Lifshitz [§]. It is 
worth noting that 

Cl^ r = v + (irn-Ho)v (2.9) 



where 



'heavy = <P K ~ V 
(CT7?) 2 \ 1 +7q 



1 + 



57777 



Q 



(2.10) 
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and H.q is a non-relativistic Hamiltonian, 



Hq= 2^ (??2 + ^ + 8^ (~ ( ^ + x ^ " ^ x ^}) + (2.11) 



well-known (in the Abelian case) from |T|, [2|]. The first term in the l/frig part is 
called the Darwin term and the second one is the convection current (spin-orbital) 
interaction. Equation ( p. 10 ) is the Foldy-Wouthuysen transformation which is nec- 



essary to keep the term linear in 7r in its canonic form. For unclear reasons the 
Foldy-Wouthuysen transformation in the context of the heavy quark theory is some- 
times referred to as "casting the Lagrangian in the NRQCD form" . For a dedicated 
discussion of the Foldy-Wouthuysen transformation in the heavy quark theory see 
Ref. |3|. 



At the one-loop level the coefficients in the effective Hamiltonian ( 2.1 1|) get some 



corrections which are specific for QCD and cannot be read off the text-book QED 
expansion. The work on calculating the one-loop coefficients for all terms through 
0(1/ nig) was recently completed [|14 |. 



The careful reader might have noticed that we consistently avoid using the term 
HQET. This is done deliberately. Although, in principle, HQET is one of a few con- 
venient technical tools for the heavy quark expansion (an alternative useful approach 
is provided by NRQCD fl5f), in many standard presentations an additional assump- 



tion is made which makes the "folklore" HQET incompatible with the OPE-based 
expansions. What is missing in the commonly accepted version is appreciation of 
the fact that all coefficients and all operators in the Lagrangian, individually, are \i 
dependent - as is quite obvious from the consideration above. The Lagrangian can 
be applied to the perturbative calculation of the Feynman graphs with heavy quarks 
where the characteristic virtual momenta flowing through all lines in the graphs are 
less than /i. The contribution of all virtual momenta above fi is explicitly included 
in the coefficients of the effective Lagrangian. 

To avoid confusion we suggest, from now on, to use distinct notations for the 
fundamental parameters of the Wilsonian Hamiltonian, on the one hand, and (per- 
turbatively defined) HQET parameters, on the other. The latter correspond to 
tending /j — > 0, after "appropriate subtraction of perturbation theory". It is clear 
that the subtraction cannot be carried out consistently in all orders, so that the 
parameters defined in this way are to be used with caution and reservations. Since 
they are widely exploited in the current literature it is convenient to introduce spe- 
cial notations. In particular, following Ref. Jl6| , we will denote the HQET version of 
the kinetic expectation value by — Ai. Formally, — Ai is given by the same Eq. ( |1.2| ) 
as fil; the perturbative contribution below \i is subtracted, however, from — Ai, 
which results in ambiguities. A pragmatically oriented reader, who is uninterested 
in the discussion of the subtraction ambiguities, to be presented below, may just 
assume that the subtraction is done, say, at one-loop order (or at two-loop order, if 
calculations are carried out at the level 0(a 2 s )). Even so, one should realize, that 
all general results obtained from the QCD equations of motion, valid for /x^ and 
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other parameters in the Wilsonian approach, are, generally speaking, inapplicable 
to the parameters of HQET. Let us parenthetically note that another key parameter 
appearing in the heavy quark theory, the chromomagnetic operator, 

/4 = 7^—(H Q \Q l -vGQ\H Q ) = --J—(H Q \QoBQ\H Q ) + 0{l/m Q ) , (2.12) 



depends on \i logarithmically, via its hybrid anomalous dimension [jl^l ; the limit 
— ► is never attempted in this case, of course. In the nomenclature of Ref. [16 
the parameter jjlq is 3A2. 

The mass formula relating uiq to the hadronic mass is a useful and, perhaps, 
the simplest application of the expansion outlined above. The 1/itiq corrections to 
the hadron mass is nothing but the expectation value of the effective Hamiltonian 



- 1 (H Q \n 2 + aB\H Q ) , , (/4 ~ Hg)h q , (0 ^ 

M Hq =m Q + A + - — + ... = m Q + A + ^ + ... 2.13 

2m %M Hq 2m c 



If we keep only the terms up to 1/rriQ, it does not matter whether the state Hq 
we average over is the asymptotic state (corresponding to mq = 00) or the actual 
physical heavy-flavor state. The difference becomes noticeable only at the level 



I/tjiq] it is expressed in terms of a few non-local correlation functions, see Ref. JlC 
for further details. 



The parameter A appearing in Eq. ( [2.13|) was introduced as a HQET constant in 



it is associated with those terms in the effective Lagrangian £(/i) (disregarded 
so far) which are entirely due to the light degrees of freedom. Needless to say that 
in the Wilsonian approach A is actually /1 dependent, A(/i). Wherever there is the 
menace of confusion the /1 dependence of A will be indicated explicitly. 

There exist a few alternative expressions for this parameter. Let us quote the 
one relating A to the gluon anomaly in the trace of the energy-momentum tensor, 

A = —^—{H \^^-G 2 \H ) mn ^ oo . (2.14) 



This expression was obtained in Ref. [TO]. Some subtleties left aside here are dis- 
cussed in detail in that paper. It is an obvious counterpart of a similar relation for 
the nucleon mass 



M N = -±-(N\^G 2 \N) . 
2M N 4a s 

The renormalization properties of the operator are quite different, however, in the 
sector of QCD with the heavy quark. For simplicity above the light quarks are taken 
as massless; introduction of the light quark masses changes only technical details. 
If the mass term of the light quarks is set equal to zero the light quark fields do not 
appear explicitly in the trace of the energy-momentum tensor. 
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3 Basic Parameters of the Heavy Quark Expan- 
sion 

3.1 The Heavy Quark Mass 

This section is a brief review of the the present status of determination of the 
heavy quark masses. An internally consistent definition of the heavy quark mass 
is of utmost importance for 1/mq expansions. Numerical reliability is essential if 
one wants to extract accurate values of the Cabibbo-Kobayashi-Maskawa (CKM) 
parameters \V C b\ and \V u b\. While these remarks are obvious (in hindsight), the 
theoretical implications were at first not fully appreciated. 

For quarks - confined degrees of freedom - there exists no natural mass defi- 
nition unlike with asymptotic states. Computational convenience and theoretical 
consistency are the only guiding principles. Convenience suggests employing the 
pole mass. However, as will be explained below, the pole mass does not allow the 
consistent inclusion of nonperturbative (power-suppressed) corrections. 

3.1.1 What, after all, is the heavy quark mass? 

In quantum field theory the object we begin our work with is the Lagrangian for- 
mulated at some high scale Mq. The mass mo is a parameter in this Lagrangian; it 
enters on the same footing as, say, the coupling constant a s with the only difference 
being that it carries dimension. As with any other coupling, mo enters in all ob- 
servable quantities in a certain universal combination with the ultraviolet cutoff Mq 
(in renormalizable theories). Although this combination is universal, its particular 
form is scheme and scale-dependent. 

The mass parameter m by itself is not observable. For calculating observable 
quantities at the scale fi <C M it is usually convenient to relate m to some mass 
parameter relevant to the scale //. For instance, in quantum electrodynamics at low 
energies (i.e. E <C m e ) there is an obvious "best" candidate: the actual mass of 
an isolated electron, m e . In the perturbative calculations it is determined as the 
position of the pole in the electron Green function (more exactly, the beginning of 
the cut). The advantages are evident: m e is gauge-invariant and experimentally 
measurable. 

The analogous parameter for heavy quarks in QCD is referred to as the pole 
quark mass, the position of the pole of the quark Green function. Like m e it is 
gauge invariant. The idea of using this parameter in the quark-gluon perturbation 
theory apparently dates back to Ref. |20|. Unlike QED, however, the quarks do not 
exist as isolated objects (there are no states with the quark quantum numbers in 
the observable spectrum, and the quark Green function beyond a given order has 
neither a pole nor a cut). Hence, m po \ e cannot be directly measured; m po i e exists 
only theoretical construction. 

In principle, there is nothing wrong with using m po \ e in perturbation theory where 
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it naturally appears in the Feynman graphs for the quark Green functions, scattering 
amplitudes and so on. It may or may not be convenient, depending on concrete goals. 

The pole mass in QCD is perturbatively infrared stable, order by order, like in 
QED. It is well-defined to every given order in perturbation theory. One cannot 
define it to all orders, however. Nonperturbatively, the pole mass is not infrared- 
stable. Intuitively this is clear: since the quarks are confined in the full theory, the 
best one can do is to define the would-be pole position with an intrinsic uncertainty 
of order ~ Aq C d • 



The issue of the appropriate quark mass has two aspects. At the perturbative 
level it is desirable to define the running mass rriQ^fi) at a scale relevant to the scale 
in the processes at hand, /i, which would incorporate contributions of all virtual 
momenta (of quarks and gluons) from M down to fi. This would allow eliminating 
unwanted large coefficients which always emerge under unsuitable choice of expan- 
sion parameters. 

If analysis is aimed at higher (power in mn 1 ) accuracy, the use of mass parameters 
that can be defined to that accuracy becomes mandatory. Whatever parameter we 
choose, it will not coincide with the mass of any observable particle, and is bound 
to remain a theoretical construction. If the construction is consistent, however, 
and allows us to achieve the desired degree of accuracy in the predictions for the 
observable quantities, that is all we need. The most straightforward choice in the 
theory with the confined quarks is as follows. We start from the original Lagrangian 
at Mo, and evolve it down to the scale /i. The original mass parameter mo evolves 
accordingly. The mass parameter appearing in the effective Lagrangian C(fi) is 
m Q (fi). 

To any order k the relation between the pole mass and the running short- distance 
one rriQ(n) takes the form Q 



(k) 

m pole 



^)Ec n W NV, cb = i. (3.1; 



Three questions (more exactly, three facets of one and the same question) natu- 
rally arise in connection with Eq. ( |3 . 1| ) : 

i) Does the series (|3.1|) converge to a reasonable number? 

ii) How well does it behave numerically in low orders, say, in the second or third 
order in a s which we usually deal with in actual calculations? 

iii) What happens when we start analyzing nonperturbative corrections sup- 
pressed by powers l/m^? 

The answer to the first question is negative, which entails negative consequences 
in the second and the third question, as we will see shortly. 

Let us begin with the divergence of the perturbative series in Eq. (|3.1| ). Al- 
though more than one reason for the factorial divergence may exist, it is easy to 

9 (k) 

In fact, Tip ole depends on \i since the a s series is truncated at the order k. This dependence 
shows up, however, only at the level O (a^ 1 ). 
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Figure 1: Perturbative diagrams leading to the IR renormalon uncertainty in rn^ Xe 
of order Aqcd- The number of bubble insertions in the gluon propagator can be 
arbitrary. 

identify a particular source, the so-called infrared renormalon singularity in the per- 
turbative expansion of the pole mass fl22| , |23H , see diagrams in Fig. 1. The bubble 
chain generates the running of the strong coupling a s . To leading order, it can be 
accounted for by inserting the running coupling constant a s (k 2 ) in the integrand 
corresponding to the one-loop expression. In the non-relativistic regime, when the 
internal momentum \k\ <C tuq, the expression is simple, 

4 r d A k A7ca s (-k 2 ) 4 rd 3 ka s (k 2 ) 
6mQ ~ S J (270***0 k 2 = 3 / A^^T- (3 ' 2) 



Now, expressing the running a s (k 2 ) in terms of a s (fi 2 ), k 2 < fi 2 

-1 

4tt " "> 2 J 3 ' 3 



a s (k 2 )= as ^ 2 ){l + ^^b\n^-\ , 6 = Hjv e _2 , (3.3) 



and expanding a s (k 2 ) in a power series in a s (/i 2 ), it is easy to find the (n + l)-th 
order contribution to Stuq, 

m Q 3 7r \ 2n J 

Observe that the coefficients grow factorially and contribute with the same sign. 
Therefore, one cannot define the sum of these contributions even using the Borel 
transformation. The best one can do is to truncate the series judiciously. An optimal 
truncation leaves us with an irreducible uncertainty ~ O(Aqcd) |23fl - 

Thus, the perturbative expansion per se anticipates the onset of the nonpertur- 
bative regime (the impossibility of locating the would-be quark pole to accuracy 
better than Aqcd)- Certainly, the concrete numerical value of the uncertainty in 
m po ie obtained through renormalons is not trustworthy. The renormalons do not 
represent the dominant component of the infrared dynamics. However, they are a 
clear indicator of the presence of the power-suppressed nonperturbative effects, or 
infrared instability of m po i e ; the very fact that there is a correction 0(Aq CD /7«q) is 
beyond any doubt. 

One can illustrate the emergence of this correction in the following transparent 
way. Consider the energy stored in the chromoelectric field in a sphere of radius 



11 



R ^> 1/rriQ around a static color source of mass tuq, 



S£coni(R) oc [ d 3 xE(i oul oc const - _ (3 5) 

Jl/m b <\x\<R 7T R 



The definition of the pole mass amounts to setting R — > 00; i.e., in evaluating the 
pole mass one undertakes to integrate the energy density associated with the color 
source over all space assuming that it has the Coulomb form. In real life the color 
interaction becomes strong at R ~ 1/Aqcd; at such distances the chromoelectric 
field has nothing to do with the Coulomb tail. Thus, one cannot include the region 
beyond Rq in a meaningful way. Its contribution which is of order Aqcd ; thus, has 
to be considered as an irreducible uncertainty which is power-suppressed relative to 

m Qi 

*n^ = fA< m \ (36) 
m Q V m Q I 

It is worth noting that the pole mass was the first example where a quantity which 
is perturbatively infrared- stable was shown not to be stable nonperturbatively at the 
level Aq CD . The observation of Refs. [2^, |23[ gave impetus to dedicated analyses of 
other perturbatively infrared- stable observables in numerous hard processes without 
OPE, in particular, in jet physics. Such nonperturbative infrared contributions, 
linear in Aqcd/Q were indeed found shortly after in thrust and many other jet 



characteristics (for a review and a representative list of references see e.g. p4|j). 

To demonstrate that the problem of divergence of the a s series for mg 16 is far 
from being academic, let us examine how the perturbative contribution to the b 
quark mass looks numerically: 

m P° le = m b (l GeV) + <5m pcrt (< 1 GeV) ~ 

4.55 GeV + 0.25 GeV + 0.22 GeV + 0.38 GeV + 1 GeV + 3.3 GeV + (3.7) 

where mfc(l GeV) is the running mass at [A — 1 GeV, and <5m pert (< 1 GeV) is the per- 
turbative series taking account of the loop momenta from 1 GeV down to zero. It is 
quite obvious that the corrections start to blow up already in rather low orders! Ex- 
pressing observable infrared- stable quantities (e.g. the inclusive semileptonic width 
B — > X u iv) in terms of the pole mass will necessarily entail large coefficients in the 
a s corrections compensating for the explosion of the coefficients in m£ ole . We will 
return to this point later (Sect. 8.1.3). 

Summarizing, the pole mass per se does not appear in OPE for infrared-stable 
quantities. Such expansions operate with the short- distance (running) mass. Any 
attempt to express the OPE-based results in terms of the pole mass creates a prob- 
lem making the Wilson coefficients ill-defined theoretically and poorly convergent 
numerically. 

Concluding this section, we make a side remark concerning the t-quark mass. 
The peculiarity of the t quark is that it has a significant width T t ~ 1 GeV due to 
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its weak decay. The perturbative position of the pole in the propagator is, thus, 
shifted into the complex plane by — §IV The finite decay width of the t quark in- 
troduces a physical infrared cutoff for the infrared QCD effects KB] . In particular, 



the observable decay characteristics do not have ambiguity associated with the un- 
certainty in the pole mass discussed above. The uncertainty cancels in any physical 
quantity that can be measured. That is not the case, however, in the position of the 
pole of the t-quark propagator in the complex plane (more exactly, its real part). 
The quark Green function is not observable there, and one would encounter the 
very same infrared problem and the same infrared renormalon. The latter does not 
depend on the absolute value of the quark mass (and whether it is real or have an 
imaginary part). Thus, in the case of top, one would observe an intrinsic uncertainty 
due to the renormalon, of several hundred MeV, in relating the peak in the physi- 
cal decay distributions to the position of the propagator singularity in the complex 
plane. 

While this fact will pose no practical problem any time soon, it will do so in the 
future, in particular when analyzing top production at linear e + e~ colliders. The 
observables can be conveniently expressed in terms of a mass m t defined at the scale 
« T t that is free from the renormalon ambiguities. It can be defined and measured 
without intrinsic limitations. It is worth trying to come to an agreement in advance 
what kind of top quark mass should be listed in the PDG tables, and thus avoid the 
problems which accompanied attempts to pinpoint the masses of c and b quarks. 

3.1.2 Short distance masses 

Since the pole mass is theoretically ill-defined, it is preferable to use a short- distance 
mass. The most popular choice is the so-called MS mass m(/x). 

The MS mass is not a parameter in the effective Lagrangian; rather it is a 
certain ad hoc combination of the parameters which is particularly convenient in 
the perturbative calculations using dimensional regular izat ion. Its relation to the 
perturbative pole mass is known to two loops j^6[ : 



< = m Q (m )|l + ^^^ + (1.566-3.73)(^) 2 + ....| (3.8) 

At \i > rriQ the MS mass coincides, roughly speaking, with the running Lagrangian 
mass seen at the scale ~ \i. However, it becomes rather meaningless at \i -C tuq. To 
elucidate the point let us note that, by definition, it adds to the bare mass m the 
same logarithmic contribution ~ a s rriQdk 2 /k 2 both at k 2 > rriQ, where it is indeed 
present, and below tuq where such a contribution is nonsensical. rfiQ(fi) logarithmi- 
cally diverges when /i — > 0. The actual contribution coming from domain below mg 
is essentially smaller, ~ a s dk 2 /k. For this reason rh(n) is not appropriate in the 
heavy quark theory where the possibility of evolving down to a low normalization 
point, /i <C rriQ, is crucial. 
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The genuine short- distance low-scale mass mg(//) at /i <C niq can be introduced 
in different ways, and it always exhibits an explicit linear // dependence, 



, — Cm 

dyU 7T 



The coefficient c m is definition-dependent. For the most natural choices \ T7[ [22| , |28 
c m varies between 4/3 and 16/9. The difference in the definition of the running mass 
in the heavy quark theory can, thus, constitute ~ 50 to 100 MeV at \i ~ 1 GeV. 
Once a particular definition is chosen, no ambiguity remains. A convenient physical 
(gauge-invariant) definition of the short- distance heavy quark mass with /i ttlq 
was outlined in |T(| (Sects. III.D and V.A). Technical aspects are discussed in detail 
in the dedicated paper [p9fl , where c m = 16/9. 



3.1.3 Which mass is to be used? 

It is legitimate to use any short- distance mass. The normalization point fi can be 
arbitrary as long as ft > Aq CD . It does not mean, however, that all masses are 
equally practical, since the perturbative series are necessarily truncated after a few 
first terms. Using an inappropriate scale makes numerical approximations bad. In 
particular, relying on mQ(m<g) in treating the low-scale observables can be awkward. 
The following toy example illustrates this point. 

Suppose, we would like to exploit QCD-like methods in solving the textbook 
problem of the positronium mass, calculating the standard Coulomb binding energy. 
To use the result written in terms of the MS mass one would, first, need to evaluate 
rh e (m e ). This immediately leads to technical problems: rh e (m e ) is known only to 
the order a 2 m e ; therefore, the "theoretical uncertainty" in m e (m e ) would generate 
the error bars ~ a 3 m e in the binding energy, i.e. at least 0.01 eV. Moreover, without 
cumbersome two-loop calculations one would not know the binding energy even at 
a few eV level - although, obviously, it is known to a much better accuracy (up to 
a 4 m e without a dedicated field-theory analysis), and the result 

M P = 2m e (0) — - (3.10) 

is obtained without actual loop calculations! 

Thus, working with the MS mass we would have to deal here with the "hand- 
made" disaster. The reason is obvious. The relevant momentum scale in this prob- 
lem is the inverse Bohr radius, hb = t'b 1 ~ am e . Whatever contributions emerge 
at much shorter distances, say, at yT l ~ "^T 1 , their sole role is to renormalize the 
low-scale parameters a and m e . If the binding energy is expressed in terms of these 
parameters, further corrections are suppressed by powers of Hb/h- Exploiting the 
proper low-energy parameters a(0), m e (0) one automatically resums chains of cor- 
rections. In essence, this is the basic idea of the effective Lagrangian approach, 
which, unfortunately, is often forgotten. 
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Needless to say, that in the processes at high energies it is the high-scale masses 
that appear directly. In particular, the inclusive width Z — > bb is sensitive to 
rrib(Mz); using MS mass normalized at \x ~ Mz is appropriate here. On the contrary, 



the inclusive semileptonic decays b — ► civ are rather low-energy in this respect p9 
and, to some extent, that is true even for b — > u. 



3.1.4 The numerical values of m c and mi, 

Accurate phenomenological determination of m& and m c at large fi (higher or com- 
parable to the quark masses themselves) is not easy since it requires very precise 
data and a high degree of control over the perturbative and nonperturbative cor- 
rections. The terms ~ (a s /7r) 2 mf, by themselves constitute typically ~ 200 MeV, 
which, thus, sets the scale of accuracy one can expect. 

The mass of the c quark at the scale ~ m c ~ 1 GeV can be obtained from the 
charmonium sum rules [Q, m c (m c ) ~ 1.25 GeV. The result to some extent is 
affected by the value of the gluon condensate. To be safe, we conservatively ascribe 
a rather large uncertainty, 

m r (m r ) = 1.25 ±0.10 GeV. 



One can argue that the precision charmonium sum rules actually fine-tune the 
charmed quark mass to a much better accuracy. The argument would lead us far 
astray, and is irrelevant for our present purposes since the convergence of the l/m c 
expansion is not good anyway. 

The desire to get rid of potentially large perturbative corrections in the b quark 
mass, scaling as m&, suggests examining the threshold region. The process e + e~ — > bb 
in the threshold domain provides us with the opportunity of determining the low- 
scale running mass, along the lines suggested in Ref. pOl. Using dispersion relations 



*<«•) = H(0) + f 2 ! (3 . u; 
1-k z J sis — q Z) 



one evaluates the polarization operator n fe (g 2 ) (and its derivatives) induced by the 
vector currents by^b, in the complex q 2 plane at an adjustable distance A from the 
threshold. Such quantities are proportional to weighted integrals over the experi- 
mental cross section; the integrals are saturated in the interval ~ A near threshold, 
and are very sensitive to the mass m&(A). 

A dedicated analysis was carried out by Voloshin j31] who considered a set of 
relatively high derivatives of ILj at q 2 = 0. On the phenomenological side they are 
expressed through moments of Rb(s), 

2tt 2 („w T f dsR b (s) 



nno) 



n\ ' • • J s n+1 



M T(S) +1) / ds Ms) exp {-(n + 1) (s - 4M 2 (1S) ) } , (3.12) 
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while the theoretical expressions for the very same moments are given in terms of 
the running quark mass and a s . The relevant momentum scale here is \x ~ m|/ \fn. 
Considering small-n moments I n one would determine rrib at the scale of the order 
rrib. The small-n moments are contaminated by the contribution of Rb above the 
open beauty threshold where experimental data are poor. Increasing n shrinks the 
interval of saturation and, thus, lowers the effective scale. On the other hand, 
we can not go to too high values of n where infrared effects (given, first, by the 
gluon condensate) explode. There is still enough room to keep the gluon condensate 
small and, simultaneously, suppress the domain above the open beauty threshold. 
In the fiducial window, /i must be large enough to ensure control over the QCD 
corrections. The latter requires a nontrivial summation of enhanced Coulomb terms 
unavoidable in the non-relativistic situation. As known from textbooks, the part 
of the perturbative corrections to the polarization operator, associated with the 
potential interaction, is governed by the parameter oi/\v \ rather than by a per se. 
In |31]] a typical scale is fi ~ 0.7 to 1.5 GeV. A sophisticated resummation technique 
was used to access the moments as high as 16 to 20. 

The crucial question is the estimate of uncertainties in determination of mj. The 
uncertainty of the fit for rrib is a meagre few MeV. This is not surprising since the 
theoretical expression for the moments very sharply depends on the value of rrib, 

I n oc e - 2n5mb/mb . (3.13) 

It is clear that the uncertainty of the fit per se is overshadowed by a systematic 
theoretical uncertainty. The calculation of the Coulomb effects in Ref. |31[] was 



not genuinely two-loop: the effect of running of a s was accounted for only through 



the Brodsky-Lepage-Mackenzie (BLM) scale fixing p2 |. Considering the fact that 
the dominant, BLM a 2 s effects were accounted for, it seems safe to assume that the 
actual accuracy of the determination of m^/i) with /i ~ 1 GeV is not worse than 
(a s (/i)/7r) 2 /i ~ 30 MeV f3~If . To be on the safe side, we increase it up to 50 MeV. This 
number includes, additionally, possible scattering in the definitions of rrib(fi) at the 
level ~ (a s /7r)/x 2 /(2mb) which can constitute 10 to 20 MeV. It is worth mentioning 
that the most labor-consuming analysis of the next-to-leading a™ +l /\v \ n effects in 



[31], although absolutely necessary for determination of a s , led to a modest change 
in the numerical value of rrib, ~ 15 MeV, as compared to the ten-year old analysis 
of Ref. HQ. 

This method could be refined in several respects to yield the overall precision in 
the ballpark of 10 MeV. At this level the accurate definition of the renormalization 
procedure is mandatory. If the experimental cross section e + e _ — > bb above the 
threshold were better measured in a reasonably wide interval, one could study lower 



3 After this review was essentially finished it was reported [Q that the analysis of the bb sum 
rules had been repeated, with some rather surprising conclusions. Resummation of the full Coulomb 
ladder is not addressed, however, in Ref. [Q. We choose to rely on the previous results until the 
matter is more closely looked at. 
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moments I n and get an accurate determination of m& without tedious analysis of 
the perturbative Coulomb corrections. 

Numerically, the b quark mass turns out to be |3"T 



m fe (lGeV) = 4.64 GeV ± 0.05 GeV. (3.14) 

The definition corresponds to the scheme in which c m = 16/9. In other words, the 
"one-loop" pole mass is 

mjp = 4.64GeV + — — • 1 GeV ~ 4.83GeV (3.15) 

9 71 

at a s = 0.336. The result is, of course, sensitive to the choice of \i. The sensitivity 
to the uncertainty in a s is around 10 MeV. The corresponding value of A(l GeV) ~ 
0.6 GeV. 

The heavy quark masses can be measured, in principle, by studying the distribu- 
tions in the semileptonic B decays [j3"5 , |36|j . Such analyses were undertaken recently 



37| , Unfortunately, the data are not good enough yet to yield a competitive de- 
termination. On the theoretical side, there are potential problems with higher-order 
corrections due to not too large energy release mb — m c ~ 3.5 GeV and/or relatively 
small mass of the c quark. In particular, the effect of higher-order power corrections 
can be noticeable. The result of analysis reported in Refs. Q is compatible 
with the value ( ^.14[ ), within the quoted uncertainties. This approach seems to be 
most natural. 

In many applications one needs to know the difference between m& and m c . If 
both masses are normalized at the same low-scale point below both masses, mb(pi) — 
m c (p) depends on p, weakly. This difference is well constrained in the heavy quark 
expansion. For example, 

rn b -m c = M B -M D + £(^-^) + ?^(±-±) + O (± 

\2m c zrrib/ 4 \mi mil \m 6 



M B + 3M B * M D + 3M D , 3 _ 3 3 

M B = ^ , M D = , p 6 = + P % 



(3.16) 



where p^ is the asymptotic expectation value of the kinetic operator, p 3 D is the 
expectation value of the Darwin term given by the local four-fermion operator and 
Pnn, Ps are positive non-local correlators HlQ| . All quantities in Eq. ( |3 . 1 6|) , except 



the meson masses, depend on the normalization point p which can be arbitrary. In 
this way we arrive at 

u 2 _ n 5 G e v 2 

m b - m c ~ 3.50 GeV + 40 MeV ^ — — — + AM 2 , I AM 2 | £ 0.015 GeV . 

0.1 GeV^ 

(3.17) 

The estimate at /i 2 = 0.5 GeV 2 appears to be in good agreement with the separate 
determinations of m& and m c from the sum rules in charmonia and T. 
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The main uncertainty in nif, — m c is due to that in the value of /x 2 . The Darwin 
term can be reasonably estimated relying on factorization ||39|| ; it is of order 0.1 GeV 3 . 
The non-local positive matrix elements and p| are expected, generally speaking, 
to be of the same order. Altogether, assuming \p i D — p 3 \ ^0.1 GeV 3 , one arrives at 
the uncertainty in rrib — m c due to the higher-order terms AM2 quoted above. 



3.2 The chromomagnetic and kinetic operators 

The non-relativistic Hamiltonian of the heavy quark ( p.ll|) in the order 1 / itlq con- 
tains two operators ( |1.2| ) and ( |2.12| ). Their expectation values in the heavy meson 
B are the key players in many applications. Let us note that the expectation values 
of the chromomagnetic operator in B and B* are related, 

/4 = l^- B (B\b l -a, u G» u b\B) ^ -^-(B^-a^blB*) . (3.18) 



The value of \Xq is known: since 



2m Q ^2 

describes the interaction of the heavy quark spin with the light cloud and causes 
the hyperfme splitting between B and B*, it is easy to see that 

H% ~ - 2m b (M B , - M B ) ~ - {M%, — M%) « 0.36 GeV 2 . (3.19) 

In the above relations the operators are normalized at the scale ~ m&. Evolving 
them perturbatively to the normalization scale /x ~ 1 GeV slightly enhances the 
chromomagnetic operator and the value of /^(/i), but this effect is numerically 
insignificant, and we will neglect it. 

In Sect. 3.1.4 we mentioned that the heavy quark mass can be extracted from 
the data on semileptonic B decays. By the same token, the value of the kinetic 
operator is, in principle, directly measurable in experiment. However, it remains 
rather uncertain at the moment. First attempts at extracting it from the semilep- 
tonic distributions were reported |37|, [38]]. Although the method seems to be the 



most promising, so far the outcome is inconclusive. Typically, the results fall in the 
interval 0.2 to 0.3 GeV 2 . We hasten to add, however, that the difference between /x 2 
and — Ai was not fully untangled. The normalization point /x should be explicitly 
introduced in the perturbative corrections to reveal the distinction between these 
two parameters. Until this is done it is unclear to which particular definition of 
— Ai the number quoted refers. As we will see shortly, numerically the difference 
between /x 2 and — Ai is of order 0.1 GeV 2 (— Ai is smaller than /x^). Note also that 

4 This remark applies also to other numerical determinations discussed below in this section, 



with the exception of Eq. (3.21) 
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the numerical analysis of the higher-order power corrections [EO] led the authors to 
the conclusion that these effects are significant in the application to the semileptonic 
spectra. 

Historically, the first attempt to determine the average of the kinetic operator 
from the QCD sum rules was Ref. HI] where a negative value ~ —1 GeV 2 was 



obtained! (This result was later revised by the author, 
thoughtful application of the QCD sum rules yielded 
prediction was later refined by the authors 



Shortly after, a more 
0.6 GeV 2 . The 



(0.5 ±0.15) GeV 2 



(3.20) 



Meanwhile a model-independent lower bound was established |!| ^3], |TU| 

fil > A G ^0AGeV 2 (3.21) 

which constrained possible values of /z 2 . It is worth emphasizing that the inequality 
takes place for /z 2 normalized at any point fi, provided fiQ is normalized at the same 
point. For large /i it becomes uninformative; so, it is in our best interests to use it 
at fi = several units xAq CD ^ 1 GeV. The inequality does not necessarily hold for 
— Ai. Exact QCD inequalities of this type will be discussed in the next section. 

Recently, an approach combining the QCD sum rules and the virial theorem was 
exploited ^6f for determination of /i 2 . The corresponding analysis is claimed to 
produce a surprisingly precise value 



0.1 ±0.05 GeV 2 



(3.22) 



noticeably below the range obtained in other analysis. We will dwell on various 
aspects of this approach, in an attempt to reveal reasons explaining the discrepancy, 
in the Sect. 6. 

The expectation value of the kinetic operator was also estimated in the quark 
models, with a spectrum of predictions: the relativistic quark model [47] gave 
about 0.6 GeV 2 or even slightly larger; the estimates of Ref. |IS| yield a close value 
0.5 GeV 2 . Attempts of extracting /z 2 on the lattices are also reported ^9 . 



The scattering of the numerical results above should not surprise the reader too 
much. As was mentioned above, although the formal definition of /i 2 and — Ai is 
the same, see Eq. (|i~2T ), these quantities are not identical, and one should expect 
noticeable differences. What is surprising is the discrepancy obtained under one 
and the same definition and within basically the same method. The most notable 
example of this type is the sum rule rule analysis, Balled al. [|2] versus Neubert p6 |. 
We will dwell on this point in Sect. 6.2. Now we would like to make several remarks 
regarding the difference between /z 2 and —\\. Conceptually the problem with Ai is 
similar to the problem encountered in the attempt to define the pole mass (more 
exactly, A H qet = M B - (m b ) pole ). 

If /x 2 is defined in the context of the Wilsonian OPE, and incorporates all soft 
fluctuations, including perturbative fluctuations from to fi, HQET assumes the 
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existence of a "genuinely nonperturbative" parameter 

- A: = /4(/i) - "(/4(//)) pert " (3.23) 

where (/z 2 (/z)) pert is a "perturbative part" of /z 2 coming from the domain below ji. As 
was emphasized long ago |50| , the procedure of separating out "perturbative parts" 



of the condensates cannot be carried out in general, simply because the notion of 
perturbation theory below fi is non-existent. The best one can do is to construct 
— Ai, starting from /i 2 (/i), in the given (say, the first or the second) order in a s (/i). 

The Wilsonian kinetic average /z 2 is fi dependent. For example, to the first order 
in a s one has 

A 4^=c,^ (3.24) 

where under some natural choice of the cut-off (i.e. the normalization point) c n = 
4/3 [0. 

Equation ( |3.24|) implies that to the first order in a s 



-A? } = Cw % 2 (3.25) 

7T 

where a s is the coupling constant (which does not run in one-loop calculations) and 
the superscript (1) indicates the one-loop nature of Ai introduced in this way. 

The first observation is obvious: — Aj is smaller than /i 2 (/i). How much smaller? 
The answer depends on the choice of /i; it is also affected by the value of a s in a par- 
ticular one-loop calculation. If fi = 0.5 GeV the values of the kinetic operator quoted 
in [[42], |46[ must be increased by approximately 0.05 GeV 2 to translate them from 



— Ai to /i 2 ; at ji = 0.7 GeV the shift constitutes ~ 0.1 GeV 2 . A similar (probably, 
somewhat larger) adjustment applies to empiric determinations from the semilep- 
tonic spectra |37j [38|. These shifts lie within the error bars of the corresponding 
analyses. 

Constructing a general "ready-to-use" cut off procedure for the Wilsonian oper- 
ators normalized at /i is a non-trivial problem going beyond the scope of this review. 
A brief discussion in heavy quark theory can be found in El], |52| and on lattices in 



Ref. [p3| ; the issue will be further elaborated in a dedicated publication j54| . Here 
we mention a few basic points relevant to the issue of /i 2 . 

A consistent field-theoretic definition of the operator Q(iD) 2 Q was given using 



the small velocity (SV) sum rules |Tfl, |5"If. Limiting oneself to the purely per- 



turbative level, one can readily outline other cut-off schemes allowing to calculate 
the perturbative /i dependence of /i 2 (/i). In doing so one must keep in mind that 
some purely perturbative cut-off schemes that have no transparent physical mean- 
ing violate the general quantum-mechanical properties of the operator in question. 
Therefore, one must exercise extreme caution in connecting various results regard- 
ing the condensates obtained in different approaches with each other. Moreover, 
perturbative cutoff procedures do not define the operators at the nonperturbative 
level. 
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Keeping in mind all these problems, one may ask whether the controversy in the 
current literature is to be taken seriously at all. The numerical comparison given 
above seems to show that the existing scattering in the value of the kinetic operator 
goes beyond the mismatches of different definitions. A better understanding of 
theoretical uncertainties in the sum rule and empiric determinations is badly needed. 
This task obviously must be the subject of a special investigation. In the absence of 
such an investigation we will focus in this review on a class of results which seems 
to be theoretically clean: exact QCD inequalities. 



4 Exact Inequalities of the Heavy Quark Theory 
in the Limit trq — > oo 

Exact inequalities reflecting the most general features of QCD (e.g. the vector-like 
nature of the quark-gluon interaction) have been with us since the early eighties 
|55fl . The advent of the heavy quark theory paved the way to a totally new class 
of inequalities among the fundamental parameters. As with the old ones, they are 
based on the equations of motion of QCD and certain positivity properties. All 
technical details of the derivation are different, however, as well as the sphere of 
applications. 



The most well-known example is Eq. ( |3.21| ) which was subject to intense scrutiny. 



As a matter of fact, this expression is just one representative of a large class, in- 
cluding, among others, the Bjorken and Voloshin relations, and the so-called third 
(BGSUV) sum rule. Since the corresponding derivations are scattered in the litera- 
ture, it makes sense to give here a coherent presentation of the subject. 
The starting point for all inequalities is the set of the sum rules 

^)-i = EI^I 2 + 2EI^I 2 . 

^ n m 

A(f,) = 2(Y,e n \T^ + 2Y,e m \T$A , (4.2) 

\ n m / 

^ = E4l^| 2 + 2Edk 3 ( / m 2 ) l 2 ) (4-3) 

" n m 

^) = -2E4K ( ;il 2 + 2E4k 3 ( ; ) | 2 , (4.4) 

" n m 

^ = E4l^l 2 + 2Edl< 2 } l 2 , (4-5) 

" n m 

- ^# = -2E4.|r 1 ( / 1r + 2E4li m 2 ) | 2 , (4-6) 
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a sequence which can be readily continued further. Here is the excitation energy 
of the k-th intermediate state ("P-wave states" in the quark- model language), 



e fc = M H {k) - M Pq , 

while the functions Ty\ and t^? 2 describe the transition amplitudes of the ground 



state B meson to these intermediate states. We follow the notations of Ref. [pc 

1 



(H% /2) \A,\P q ) = -t 1/2 ( Vi -v 2 ),, (4.7) 



and 



{H% 2) \A,\P Q ) = -J-ir^e^^V! + v 2 f{ Vl - v 2 y , (4i 



where 1/2 and 3/2 mark the quantum numbers of the light cloud in the intermediate 
states, j n = l/2 + and 3/2 + , respectively, and is the axial current. Furthermore, 
the slope parameter p 2 of the Isgur-Wise function is defined as 

2Kr( p QW\ &°Q)„ \ P Q) = 1 - £V)y + 0{f) . (4.9) 

All quantities on the left-hand side of Eqs. ( |4.1| ) - ( |4.6|) depend on the normalization 
point /i, since they are defined in the Wilsonian sense. The right-hand side also 
depends on this point. Thus, all sums over the excited states appearing in Eqs. 
( (4.1[ ) - (|4.6p imply a cut-off at the excitation energy equal to /i, 

e k <fi. (4.10) 

Equation (|4.1|) is nothing but the Bjorken sum rule \fi% . It was discussed in 



great detail in Ref. |56[ where, among other things, the cut-off was introduced in the 



form flj.lOf) almost a decade ago! Equation (|4.2|) was obtained by Voloshin []58| . The 



expression for ^ is the BGSUV sum rule The next one is new. The last two sum 
rules are obtained along the same lines; the one for the Darwin term p 3 D was presented 
in [60 1 . A unified field-theoretic derivation of all these relations was discussed in 
Ref. [0] where the transition amplitude in the SV limit was considered (see also 
Sect. 8.2.1). Since the presentation in that paper is sufficiently pedagogical, the 
interested reader is referred to it directly. In this section we will give an alternative 
(new) derivation of Eqs. ( |4.3|) -( |4~6l) . A new derivation of Eq. ( |4.2f ) will be sketched 
in Sect. 6.1. 

At first, however, we pause to make a few explanatory remarks regarding \x de- 
pendence. To avoid inconsistencies it must be done in the same way in all theoretical 
calculations, including perturbation theory, the definition of the heavy quark current 



5 Higher moments than those presented above are not very useful in practice since they are 
saturated at higher e, the /i dependence becomes too essential, etc. 
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(and, hence, the Isgur-Wise function) and so on. This cut-off scheme does not coin- 
cide with the one implied by dimensional regularization, although they are related in 
the perturbative domain /i 3> Aq CD . In particular, the logarithmic scale dependence 
of £> 2 (/u) is the same. The bound q 2 > 1/4 |57j holds only for this definition of the 
renormalized Isgur-Wise function and is not necessarily true for other definitions. 
Without the cut-off at = fi the sums over the excited states in the above relations 
diverge at high excitation energies. The high-energy tails of the sums are dual to the 
perturbative transition probabilities which, thus, determine the [i dependence of the 
operators at hand. Cutting the integrals over the excitation energies (i.e. limiting 
the sums by (e/ c ) max = //) is a physical way of introducing the normalization point 
without violating the general properties of QCD and not endangering the quantum- 
mechanical aspects of the field theory. It is often convenient to use the exponential 
weight e~ £ ^ in the sums. As long as the weight function is universal, all relations 
are the same as with a step-like cut-off. 

Now, we proceed to a consideration which illustrates the physical meaning of the 
sum rules presented above. 

Since in the heavy quark limit, mq ^> Aq CD , QQ pairs are not produced, we will 
use the quantum-mechanical language with respect to Q (but not the light cloud, of 
course). Moreover, it is convenient, at the first stage to assume Q to be spinless. The 
Q spin effects are trivially included later. Then the lowest-lying states, the S'-wave 
configurations corresponding to B and B*, are spin-1/2 fermions, with two spin 
orientations of the light cloud. We shall denote them \Qq); the spinor wavefunction 
of this state is ^o- It is obvious that 

2 2 

(n |Q(*£j)(i£*)Q|flo) = ^S jk *l* - ^*S<7 ifc tt = £< n ok>)(n|7r fc |n ), 

n 

(4.11) 

where a complete set of intermediate states is inserted. They are spin-1/2 states (of 
the opposite parity with respect to |fi )), generically denoted by <jM\ and spin-3/2 
states x ■ We will use the Rarita-Schwinger wavefunctions for the latter, i.e. a 
set of three spinors xi obeying the constraint uiXi = 0. The normalization of these 
spinors is fixed by the sum over polarizations A 

£Xi(A)xj(A)=<%-^ . (4.12) 

Defining the reduced matrix elements a n and b m as 

(0("V^o> = an0 {n)t ^o , (x (m) Kil^o> = 6mxf )f *o , (4-13) 
where (j)^ and x^ stand for the states as well as for their wavefunctions, we get 

Mg = -6£|a n | 2 + 2£|6 m | 2 , (4.14) 
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and 



3$> n | 2 + 2$> r , 



(4.15) 



These expressions can be immediately generalized to the actual case of the spin- 
1/2 quarks Q. The quantities a n and b m are to be understood as the matrix elements 
of biDb between the B meson and higher even-parity states. They are related to 



r, ?2 and t-To as follows: 



'1/2 



'3/2 



T 



(n) 
1/2 



and, therefore, 



On 



T. 



(m) 
3/2 



_1_6, 

V3e r , 



(4.16) 



/4 



W|2 



X! e nl T l/2 



2 I M|2 
ml '3/2 I 



/4 



(n)|2 , 9 2 |_M|2 
1/2 I e ml r 3/2 I 



Relations ([4. 16|) are most easily obtained from the fact that the SV amplitudes 
of transitions to the corresponding states are given by the overlap {n(v)\B(v = 0)) 
and by exploiting Eq. (|6.7f) for \B(v)) from Sect. 6.1. 

The sum rules ( f4.1| )-( p~6| ) obviously entail a set of exact QCD inequalities. The 
first one is nothing but the Bjorken inequality g 2 > 1/4 [[57]] which was already 
mentioned above. Others from this sequence are: 



A(^)>2A 1 L 2 ( / x) 



(4.17) 



/40) > -AxAOi) and /4(/i) > 3Af q\h) - - 



/4 = /4 + 9 



W|2 



C n l T l/2 



> Pg > 



(4.18) 
(4.19) 



p^) > Ai/i^/i) and p 3 D (p) > -A 2 A(/i) and p^(p) > 3A? ( £? 2 (/i) 



r (n)\2 
r l/2 I 



> -PlsQ*) > and p^(p) > 



(4.20) 
(4.21) 



Here Ai is the first excitation energy. Strictly speaking, the lowest excitation of B is 
Bit; the pions are coupled rather weakly, however. Moreover, the Bw intermediate 
state becomes irrelevant in the limit N c — > oo. Therefore, for all practical purposes 
we can equate Ai to the mass difference of B** and B, i.e. Ai ~ 0.5 GeV. 

It is important that the above inequalities are based only on general QCD rela- 
tions and the heavy-quark limit: no recourse has been made to model assumptions. 
In particular, dynamics of the the light cloud is included in full and exactly. The 
inequality p z D — (1/3) p\ s > following from Eqs. (|4.5f) and (|4.6|) is useful for the 
zero recoil sum rule discussed in Sect. 8.2.3. 
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The values of T3/2 and, in particular, T1/2 are not well known at present. For the 
lowest states with e ~ 0.5 GeV they were evaluated in QCD sum rules. According 



to [5T|, rj* 2 m 0.25; the similar estimate for t9X is r^X — 0.4 ±0.1 |B2|. A model 



estimate []5£| falls quite close, rfX ~ t^X ~ 0.31. Recently, the value r^X ~ 0.35 was 



derived |63|] from the experimental information on the yield of the excited charm 



states in semileptonic B decays, yet also with sizeable uncertainty. Substituting 



T i/2 ~ 0-25 in Eq. ( 4.1iJ| ) we arrive at a stronger bound on /i 2 , quoted in Eq. (4.22) 



'1/2 

below. 



All these inequalities are saturated if it is only the first excitation that contributes 



to the sums. This is explicitly so in the so-called two-doublet model of Ref. p4"| . If 
jj, is not too large, in the ballpark 0.7 to 1 GeV, the actual situation seems to be not 
too far from this model. (In the two-doublet model the higher excited states are 
exactly dual to the perturbative corrections). 

Assigning to the estimate of t^X above an uncertainty ±30%, one obtains from 

Eq. flcrgp 

/4(0.7 GeV) = (0.6 ± 0.15) GeV 2 . (4.22) 
Moreover, assuming that f\ g 2 (p) = 0.8 ± 0.2 and Ai = 0.5 GeV we get 

A(/i) = 0.35 to 0.75 (0.5 to l)GeV, (4.23) 

fxl(fx) =0.25 to 0.56 (0.38 to 0.8) GeV 2 , (4.24) 

p 3 M = 0.13 to 0.28 (0.2 to 0.4) GeV 3 . (4.25) 

The first interval assumes complete saturation by the first excitation, while the 
interval indicated in parenthesis assumes 70% saturation. 

The exact QCD inequalities presented above take place at any fj,, provided that 
all parameters refer to the Wilsonian formulation. At large \i they become non- 
informative since all dimensionful parameters (except fiQ, p\ s ) are dominated by 
the large perturbative pieces scaling like the corresponding power of /i. 

Concluding this section, let us stress again that the information encoded in the 
exact inequalities is model- independent. It makes sense to impose the corresponding 
constraints in any analysis aimed at extracting the set of the fundamental parame- 
ters, say, from the semileptonic spectra. If a particular determination leads to results 
incompatible with these inequalities, one may be sure that something is wrong either 
with data or with the analysis. 

In the practical applications discussed in the later sections we will use, as a 
central value, /z 2 = 0.5 GeV 2 normalized at 0.5 to 0.7 GeV. 

6 The experimentally measured slope is centered around 0.9, with the errors ~ ±0.2. To get 
g 2 (p) from this slope one must eliminate l/m c and hard perturbative corrections. According to 
estimates in pq, pq], this amounts to subtracting ~ 0.1 from the experimental number. 
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5 More on fi^ > {1q 

In view of importance of this particular inequality a few additional comments are 
in order. Unlike quantum-mechanical potential problems, in QCD the heavy quark 
kinetic operator is expressed in terms of the covariant derivatives i\j = —iDj = 
—idj — Aj which contain the vector potential Aj. The presence of Aj leads to non- 
commutativity of different spatial components of the momentum operator in the 
presence of the chromomagnetic field B, 

[TTj^k} = iG jk = -iejkiBi . (5.1) 

This non-commutativity immediately leads to the lower bound on the expectation 
value of 7? 2 ||39||, in full analogy with the uncertainty principle in quantum mechanics. 



The simplest quantum- mechanical formulation was suggested in Ref. f45| , 

(B\(an) 2 \B) = (B\n 2 + SB\B) > . (5.2) 

This inequality has a transparent physical interpretation. We deal here with 
the Landau precession of a colored, i.e. "charged" particle in the (chromo) magnetic 
field. Hence, one has (p 2 ) > \B\. Literally speaking, in the B* meson the quantum- 
mechanical expectation value of the chromomagnetic field is suppressed, (B z ) = 
— Hg/3- It completely vanishes in the B meson. However, the essentially non- 
classical nature of B (e.g. (B 2 ) > 3{B) 2 ), in turn, enhances the bound which then 
takes the same form as in the external classical field. 



The SV sum rules [|4], |10[ translate the quantum-mechanical derivation of the 
above inequality to the field-theoretical language; the difference /i^(/u) — is 
represented as the imaginary part of the transition operator which is given by a sum 
of certain transition probabilities. 



6 Virial Theorem 

In this section we continue the discussion of the kinetic operator from a different 
perspective. As was mentioned, the expectation value of 7f 2 in the B meson was 
estimated in the QCD sum rules directly and [HH] via the virial theorem. A 
strong discrepancy between the two determinations, which do not overlap within 
their respective error bars, prompts us to revisit the issue. Two questions arise. A 
practical one is what result is more trustworthy and why. A more theoretical prob- 
lem is what went wrong with the evaluation of the uncertainties in the competing 
calculations? 

Clearly, only the authors themselves could answer these questions in full. We 
would like to suggest, however, some qualitative insights which, hopefully, may shed 
light on the discrepancy and be helpful in future refined treatments of the issue. We 
must start, however, with the brief description of the virial relation itself. 
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6.1 Alternative derivations 



The kinetic and chromomagnetic matrix elements are defined via the expectation 
values of the corresponding operators in B meson. In Ref. |J a certain relation for 
the matrix elements of Qtt 2 Q was obtained involving a new, chromoelectric operator 
QgEQ. The transition matrix element of QgEQ between the states with different 
velocities was considered, and it was found that 

(H Q W\QiGTQ\H Q (v)) = 



2M Hq 



i (t/V* - «V") \j^-(H Q {v)\Q^Q\H Q {v)) + 0{vv - 1) j . (6.1) 

This equation is valid for the spinless Hq\ for non- vanishing spins it is valid after 
averaging over polarizations of Hq. In the rest frame of the initial hadron (v = 0) 
the relation is nontrivial for {//, u} = {0,k}; then Gok = —E%t a where E is the 
chromoelectric field. Then Eq. (|6.1|) takes the form 



1 -{H Q {v')\Q{iE*)Q\H Q m = V \\ + ... (6.2) 



2M Hq n vv ^> vv " 3 

where the ellipses denote 0(v' 2 ) terms which will be systematically omitted. 

Below several complementary pedagogical derivations will be presented. First, 
however, a few words as to why Eq. (|6.1|) is sometimes referred to as a virial theorem 
M. Assume that, instead of QCD with its messy dynamics of infinite number of 
degrees of freedom residing in the light cloud, we consider a non-relativistic bound- 
state problem: a heavy quark interacting with a lighter one via a potential. Using the 
formalism of non-relativistic quantum mechanics it is easy to rewrite the left-hand 
side of Eq. (]6.2|) as follows: 

{e- im ^' s igE\x)) -> m Q v' i {gx i E i (x)) = mQv' i {gx i ViV(x)), (6.3) 

where the angle brackets denote averaging over the given bound state (in the rest 
frame), V is the binding potential and, we restore the coupling constant contained 
in G^ v . Equation ( |6.3|) now takes the text-book form of the virial theorem, 



(ttD = m(x k d k V(x) .) (6.4) 

No sum over k is implied here. (Note that in quantum mechanics (mxkdkV(x)) 
remains the same being calculated for the heavy "quark" Q or the lighter one q.) 

Equation ( |6.1| ) might be helpful in various applications. It was exploited, in 
particular, in studying the properties of the kinetic energy operator under mixing 

The most elementary derivation of Eq. (|6.1| ) is as follows. If v — and I^T'I 1 
then 

(H Q (v')\Q(iron ~ 7T t 7r )Q\H Q (v)) = (H q {v')\Q{tt - vn? + ^)^Q\H Q {v)), 
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where all 7r's act to the right, and we have used the equation of motion neglecting 
the terms suppressed by 1/rriQ. Moreover, to the first order in v' the operators tt 
on the right-hand side can be considered as acting to the left (the difference is the 
full derivative which is obviously (9(|ir'| 2 )). When ttq — ttv' acts to the left it again 
produces zero due to the equation of motion and we are left with 

(H Q {v')\Q{itv')nQ\H Q {v)) (6.5) 

which is equivalent to Eq. ( |6.1| ). 

It is instructive to give a different derivation, more directly related to the quantum- 
mechanical aspect of the problem. Consider 

(# O (0)| W |#o(0)) = (jy O (0) Kj-TTo 1 TToTTfcl jy Q (o)> = 

(^(0)1^0 1 (ttott, - 7r fe 7r )|#Q(0)> . (6.6) 

The operator , Kq 1 {v ; k) acting on \Hq) is nothing but the generator of the boost along 
direction v: 

\H Q {v)) = \H Q (0)) + 7r l vn\H Q (0)) + 0(v 2 ) . (6.7) 

This is a useful relation nicely elucidating the meaning of the small velocity sum 
rules. It represents the first-order perturbation theory in 5H = vn; the unperturbed 
Hamiltonian is TCq = 7r (further details can be found in [0, Eq. (178) and Sect. VI). 
In the second-quantized notations the very same relation takes the form 

\H Q (v)) = \H Q (Q)) + J d 3 xQn^ l v7TQ(x)\H Q (0)) + 0(v 2 ) . (6.8) 

Keeping in mind that (HQ(0)\iTi\HQ(0)) vanishes and invoking the definition 

[ttV*] = l G ok = -iE\ (6.9) 

we immediately rewrite Eq. (|6.6| ) in the form 

(#Q(0)|g^7f7r fc Q|# Q (0)> = (H Q (v)\QiG ok Q\H Q (0)) (6.10) 

which is again the virial relation. 

The simple derivations above are formulated in terms of the on-mass-shell states. 
In some applications (e.g. the QCD sum rules), rather than dealing with the on- 
mass-shell states, one works with the Green functions and amplitudes induced by 
interpolating currents. For future applications we find it useful to formulate yet an- 
other derivation of the virial theorem based on the technique of the Green functions. 
In this approach we consider the three-point functions induced by appropriately 
chosen currents, and then use the reduction formula. The relevant correlators are 
treated in the background field method (for a review see [p9~|), although actually we 
use very little of the background field formalism. The advantage of this approach is 
its manifest field-theoretic nature. One operates with full QCD and takes the limit 
niQ — > oo, as in Ref. 0. 
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Figure 2: Three-point function relevant to the derivation of the virial theorem in 
QCD. 



Consider the three-point function depicted in Fig. 2. The central vertex in this 
triangle is the operator G° l . Two other vertices are generated by the interpolating 
currents 

J = Qi^q and = qi'j^Q , (6-H) 

where the pseudoscalar current is chosen merely for defmiteness. The sides of the 
triangle are Green's functions of the quarks in the background gluon field. The 
reduction theorem tells us that to get the transition amplitude (Hq\QG 0i Q\Hq) from 
this three-point function we amputate it: multiply by (p 2 — Mjj Q ) and (p' 2 — Mjj Q ) 
and tend p 2 and p' 2 to Mjj Q . This singles out the double pole whose residue is 
proportional to the meson-to- meson transition on G° l . 

The expression for the three-point function of Fig. 2 has the form 

*Tr Ul,-—, j G 0i — - 1 ^75^—} , (6-12) 

where tt is the momentum operator, Tr implies taking trace in the color, Lorentz and 
(abstract) momentum spaces, and p and p' are external momenta flowing through 
the currents. In the non-relativistic limit, keeping only terms linear in v', one can 
readily rewrite Eq. ( |6 .12 ) as follows: 



iTr 



G 



Oi 



(e + 7T ) - nv' (e + tt ) 



. 70 + 1-^7 1 7o + l . 1 , . 

^75 t> -i 7T~ *7s \ , 6.13 

2 jt — m q 2 J 



where e + ttlq is the time component of the external momentum p (see Fig. 2), and 
we set e' = e, since only the double pole of the type 

' (6.14) 



(e-A) 



is of interest for singling out the ground state Hq from the tower of states created 
by the interpolating currents J. 
Now, let us use the fact that 

G 0i = -t[7c°, tt'} = -i {(e + 7r o y - 7r*(e + vr )} . (6.15) 
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Substituting Eq. ( |6.15| ) in Eq. ( |6.13| ) and expanding the denominator in v' we 
observe that the square bracket reduces to 



((e + iro)** - 7r*(e + 7T )) - r — 

\ / If -+- ' 



(e + 7T )-7lV ,x - J [e + 7T 



terms with single pole+ 



1 -thTV . 1 . +0(v' 2 ). (6.16) 



(e + tto) (e + 7t ) 

In the limit e — > A, the residue of the double pole (in the second line in Eq. ( |6.16| )) 
is just 

which proves Eq. (|6.2|) . Note that it is crucial to single out the double-pole contri- 
bution. 

By the same token, a similar consideration would lead us to Eq. (94) from Ref. 
ill , and, eventually, to the sum rule (|4.2|) . Consider, instead of the three-point 



function ( |6.12| ), the two-point function 

zTr I ^ 5 in t ^ 5 "J— — r ' ( 6 - 17 ) 

{ p+ fi-rriQ f -m q ) 

with the external momentum p = itlqv^ + e M . It is assumed that t> M = {t>o, v} where 
\v\ 7^ but \v\ -C 1 and e M = {e, 0}. We then take the limit wiq — > oo, expand in 
v keeping the terms quadratic in this parameter, isolate the double pole Q6.14j) and 
compare the result with the phenomenological representation for the same two-point 
function. Equation (94) from Ref. [TO] is immediately reproduced. 



6.2 Virial theorem and QCD sum rules 

In this section we dwell on the problem of determining the kinetic energy from 
different QCD sum rules. 

The technical reason of the existing discrepancy is quite evident: in the sum 
rule analyzed by Ball et al. the free-quark diagram gives a contribution in the 
theoretical side of the sum rule. This contribution is then "corrected" by 0(a s ) 
terms and nonperturbative condensates, a standard situation in the QCD sum rules 
technology. On the other hand, Neubert suggests ||6| considering the virial partner 



of 7r 2 , and analyzing the QCD sum rule for QEQ rather than for Qtt 2 Q. In this case 
the free-quark contribution on the right-hand side is absent. The theoretical side of 
the sum rules is solely due to interaction terms which must be kept relatively small, 
by necessity. Although such a regime is much less studied than the standard one, 
some experience is still available. 

To get an idea of how the QCD sum rules work in such problems it is instructive 
to consider a toy model which was used for this purpose previously more than once 
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fTOfl . We mean illustrating the sum- rule approach in the three-dimensional harmonic 
oscillator. 

Assume that we have an infinitely heavy "quark" Q and a lighter but still non- 
relativistic "quark" q with mass m <C itlq, interacting via the harmonic oscillator 
potential, 

where r is the distance between Q and g. 

Needless to say, this model is exactly solvable. We will pretend, however, that we 
do not know it, and will calculate the matrix element of interest using the sum-rule 
approach. Parallelizing QCD we will consider the expectation value of the kinetic 
energy and its virial partner \xidiV(x). Since the potential is quadratic, the virial 
partner is the potential energy operator itself, 

2 -»2 

^ rmo r 

E pot = ^—. (6.19) 
The kinetic energy operator in the model at hand has the form 

^■■ = £ = -i A ' < 620) 

As is well-known, the expectation values of the E pot and E^ in are indeed always 
equal to each other in the harmonic oscillator; for the ground state 

Q 

(£ kin )o = (£pot)o = T • (6.21) 

Note that 

Ekm + E pot = H (6.22) 

where H is the Hamiltonian of the system. 

The important element is that, like in QCD, the potential energy vanishes for free 
particles and appears only when the interaction is included. Thus, its calculation 
is technically similar to the QCD sum rule for QEQ. The kinetic energy, clearly, is 
nonzero even for free - but moving - quarks and, in this respect, is the counterpart 
of Qn 2 Q. 

Now, in the spirit of the sum-rule approach, we pretend that the only quantities 
we are able to calculate reliably are the correlation functions at short (Euclidean) 
times. Namely, we consider the following amplitudes: at the initial moment of time 
the quarks are at zero separation, at final (Euclidean) time T = t\ + T2 they are 
at zero separation again; at t\ the insertion of the operator E^ Q or E pot is made. 
These amplitudes are perfect analogues of the Borel-transformed two- and three- 
point functions in the QCD sum rules. 

All relevant formulae are collected in Ref. |61| where a toy-model calculation in 
a related problem is carried out. The two-point function is 

S(T)=K {0 , r|0, „) = (-) j-^r,, (6.23) 
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where K(r, r|0, 0) is the amplitude (time-dependent Green function) of the propa- 
gation from the point (0, 0) to the point (r, r) in the Euclidean time, 

K(f,T\6,0) =J2r n (0)M^)e- E - T = 

n 

mUJ ( mUJ -2 w \\ In nt \ 

exp — - . — -r cosher) . (6.24) 



s 27r sinho/T / y 2sinh(u;r) 

The three-point functions with the insertion of E pot and -Ekin are 

S kin = J d 3 rK(0, n + r 2 \r, n) E kin K(r, n|0, 0) , (6.25) 

plus the same expression with the insertion of E pot . Following the standard routine 
of the QCD sum-rule practitioners we will analyze the sum rules at the symmetric 
point, 7~i = r 2 = T/2. Using Eq. ( |6.24| ) above it is easy to find that 



3 cosh(wT) + 1 3 cosh(wT) - 1 

S kin = -ujS — . S pot = -ujS — . , 6.26 

4 sinn(wi ) 4 sinn(wi ) 

where S is defined in Eq. ( |6.23| ). 

Let us review how the sum rule technology works. As a classic example, we 
determine the ground state energy of the three-dimensional oscillator. To this end 
we will analyze two sum rules 



'niT^V 2 ( uo 2 T 2 19cu 4 T 4 

n=0 

and 



^kin + ^pot = EV'n(0)C(0)Ene 
n=0 



^f^y /2 f 1+; £!i!-i^i! + ..y (6 . 28) 

2 1 2vr / 1 12 288 / K ' 



Equation ( |6.27| ) is a two-point function sum rule, Eq. ( |6.28| ) is a three-point func- 



tion sum rule with the insertion of the Hamiltonian. The right-hand sides in both 
expressions are given in the small-T expansion since we pretend not to know the 
exact solutions. Note, that the leading terms in these expansions are cu-free; they 
correspond to the propagation of free "quarks" and are analogous to the free-quark 
diangle and triangle diagrams in QCD. The subsequent terms are proportional to 
powers of to and are analogous to the condensate corrections in the sum-rule ap- 
proach, and, indeed, can be easily calculated without the exact solution at hand. 

To single out the ground state contribution it is desirable to go to as high values 
of T as possible since all higher states are exponentially suppressed at large T. 
We cannot go to the asymptotically large values of T, however, since the power 
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corrections blow up. Upon inspecting Eqs. ( 6.27|) and ( |6.28| ) it becomes clear that 



if we want to keep the power corrections under control (i.e. the last correction kept 
to be less than ~ 30%) we can not go higher than T ~ 1.5a; -1 . At such values of 
T the suppression of the higher states in the sum rules, although quite visible, is 
still not good enough. Thus, the relative contribution of the first excited state in 



Eq. (|Q7|) is (3/2)e~ 2wT , while in Eq. flOSQ its relative weight is (7/2)e -2wT . We 



can see it in a different way: taking the ratio of Eqs. ( |6.28|) to ( |6.27|) would produce 



3o>/2 if all higher excitations were fully suppressed. Instead, we get ~ 0.75(3o>/2) 
at T ~ 1.5a; -1 . 

If we aim at better accuracy we have to have at least a rough idea of the con- 
tribution coming from the excited states. The standard strategy in the sum rule 
analyses is representing all states higher than the ground state by the free-quark 
approximation, starting from some effective continuum threshold e c . The particular 
form of the corresponding spectral densities is determined by the free-quark term in 
the power expansion. In the sum rules ( |6.27|) and ( |6.28| ) we have, respectively, 



/ m \ 3 / 2 9 r°° 



e CT 



and 



EW0X(0)^ - \ (£) 3/2 ^/f^ 3/2 ^ T - (6-29) 

Transferring the "continuum" contribution to the right-hand side of the sum rules 
we arrive at 



+ ... 



480 

and 

wv ; ^ 0V ; 2 \2nJ \3^Jo 12 288 

(6'.30) 

Assume that our idea of the excited state contributions is roughly correct. If 
T can be chosen large enough so that the higher states are sufficiently suppressed, 
on the one hand, and small enough to keep the power expansion in the right-hand 
side under control, on the other, the ratio of the two expressions must be close to 
Eq = (3/2a>). The interval of T satisfying the above conditions is called window 
(fiducial domain). Of course, the constancy (T independence) of the ratio is valid 
only to the extent one can neglect the higher-order power corrections and the error 
of the continuum approximation. Thus, we expect the ratio to be approximately 
constant. The solid curve in Fig. 3 shows the ratio, as a function of T in the fiducial 
domain (the plot corresponds to e c = 2.5a;; the value of the "continuum threshold" 
can be fitted in its turn). We see a clear-cut plateau, and the height of the plateau 
differs from the exact ground state energy (Eq = 3o>/2) by at most ±5%. 
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Figure 3: Sum rules for the three-dimensional harmonic oscillator. The solid curve 
presents (H) (in units of |a>) vs. T. The dashed curve gives (E kin ) in units of |cj, 
the dotted curve {E pot ) in units of |a\ The horizontal axis is T in units of u~ 1 . The 
exact results for the ground state energy, kinetic energy and potential energy are 
|u;, jUJ, jou, respectively. 



So far the application of the sum rule technique is successful. What happens, if 
instead of the sum of the kinetic and potential energy operators we would analyze 
them separately? 

From the exact solution of the problem we know that the most drastic change 
is the occurrence of the off-diagonal transitions in the three-point sum rules for S^n 
and Spot- These off-diagonal contributions are absent, by construction, in S and are 
obviously absent in the sum rule for (H), too, since the off-diagonal matrix elements 
of the Hamiltonian vanish. On the contrary, the off-diagonal matrix elements of the 
operators E kin and E pot do not vanish and are numerically relatively large compared 
to the diagonal matrix elements. For instance, (r 2 ) i = — 2(r 2 )oo- Moreover, if 
the first excited state appearing in S or (H) has the excitation energy 2uj (above 
the ground state), the first excited state in (E^m) and (E pot ) is separated from the 
ground state by only u>, as explained in detail in Ref. |JT|. The fact that the excited 
states lie closer to the ground state and their relative weight is relatively large makes 
it harder to suppress them sufficiently/and or approximate them by the free-quark 
continuum. The fiducial domain is bound to shrink. 

A closer inspection shows that the situation is especially bad in the sum rule 
for (Ep t), since in this case the continuum simply vanishes in the free quark ap- 
proximation. This is obvious from the definition of the operator, and is clearly seen 
from the small-T expansion of the second expression in Eq. ( |6.26| ). Following the 
standard strategy we would be forced merely to neglect all higher states. The exact 
solution shows |JT] that the excited states contribution consists of two components 
- the diagonal transitions and off-diagonal ones - and each component separately 
does not vanish in the free quark approximation and is large, but they have opposite 
relative signs. For highly excited states the sign-alternating terms, after summation, 
are smeared to zero, so this approximation of "no continuum" is not bad. The first 
excitation, however, is too conspicuous and is not annihilated by the second one. For 
values of T up to ~ 1.5a;~ :L the first excitation plays an important role in the sum 
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rule, and effectively screens the ground state contribution. Accepting the "stan- 
dard" continuum model (i.e. no continuum in the case at hand) we significantly 
underestimate the expectation value of E pot in the ground state. Figure 3 illustrates 
the trend. Assuming that the fiducial domain is the same as in the sum rule for S 
or (H) we get around 0.5 of the exact value. As a matter of fact this means that 
the fiducial domain shrinks to zero in this case. 

Analysis of the kinetic energy operator is better, but not much better. The 
diagonal and off-diagonal transitions are the same in the absolute value as above 
(this is obvious from Eq. (|6.26 )), but the relative sign is now positive. The free 
quark approximation gives a non-vanishing continuum contribution. This is always 
good news. The bad news is that the lowest lying excited state (the off-diagonal 
transition), even after smearing, is not well-described by the free-quark curve. On 
average, the free-quark curve falls below the actual value of the (smeared) spectral 
density. To make a relatively precise prediction we would have to go to such high 
values of T where the power expansion, truncated at a couple of the first terms, 
fails. If we stay inside the former fiducial domain, we with necessity overestimate 
the value of the kinetic energy in the ground state. The outcome of the sum-rule 
calculation includes a positive and rather significant contamination from the first 
excitation. The extent of this contamination is clear from Fig. 3. Instead of 3a>/4 
we get roughly 1.5 • 3u;/4. The sum is very close to the exact result, 3o>/2, but since 
the divergences go in the opposite directions, the difference, instead of being zero in 
accord with the virial theorem, is quite large. 

What lessons can be inferred for the QCD calculations? Although the toy model 
discussed above is far too simple to closely follow the QCD pattern, some features 
are still rather close. The virial partner of n 2 in QCD is E. It is quite obvious that 
the sum rule for E does not have the free-quark graph. In other words, just like 
in the case of (E pot ) the free-quark continuum vanishes H4~6f . This was believed to 
be a positive feature of the analysis. We see now that it is not. We expect a large 
physical spectral density in close proximity to the ground state contribution (caused 
by off-diagonal transitions), and of the opposite sign. To suppress it one would 
need to go to such high values of the Borel parameter where - alas - we cannot go 
without losing control over the power expansion. The screening of the ground state 
contribution by the off-diagonal transitions results in an underestimated value of 
(gE), which through the virial theorem leads to an underestimated value of 7? 2 . 

At the same time, and for the same reason, the calculation of Ref. |L2]j , which 
follows the standard pattern, is expected to give an overestimated prediction for 
7? 2 . Whether the genuine value of fi 2 is closer to the one edge of the interval or to 
the other - we do not know. It seems clear, however, that the error bars in Eqs. 
( |3.2(J| ) and ( |3.22| ) reflect only uncertainties in various parameters entering the sum 
sum-rule calculations and do not reflect the fact that the continuum saturation in 
the problem at hand is likely to deviate from the standard pattern.^] 

7 The authors of Ref. estimated the uncertainty of the continuum model reflected in 
Eq. (pT20l) at the 30% level, i.e. 0.15 GeV 2 . 
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The actual value of tt 2 lies, most probably, somewhere in between the two esti- 
mates. Unfortunately, the potential interaction is strictly speaking non-existent in 
QCD, let alone the quadratic potential of the oscillator model. Therefore, mechani- 
cally averaging the estimates ( |3.20 ) and ( 3.22 ) is unlikely to suppress the off-diagonal 
transitions and yield a more reliable approximation to the genuine value of fi 2 . In 
no way do we suggest such an averaging. It is clear that more sophisticated models 
for the off-diagonal transitions are needed in both calculations. Once we realize that 
the standard continuum approximation is inadequate, developing a better one does 
not seem to be a hopeless problem. 



In summary, it seems safe to say, that the result obtained in Ref. [[4(J signifi- 
cantly underestimates the kinetic operator. The sum rule technology for tt 2 seems 
to operate in a somewhat better environment than that in the case of E. This con- 
clusion is indirectly supported by the fact that the alternative analysis relying on 
the exact QCD inequalities (plus the data on the slope of the Isgur-Wise function 
and the observed spectrum of the excited states) favors larger values of /i^.. 



7 OPE for Inclusive Weak Decays 

The number of applications of the heavy quark expansion dramatically grows when 
one includes external (non-QCD) interactions of the heavy quark, e.g. electromag- 
netic, weak and so on. A variety of problems can be formulated in the language of 
the OPE - inclusive heavy quark decays are the best-studied example of this kind. 
In this case, the operator product expansion is built in a slightly different way. The 
expansion parameter is not necessarily l/mg, but rather it is regulated by the en- 
ergy release in the problem at hand, or by other external parameters. The general 
consideration runs parallel to the treatment of <r(e + e~ — > hadrons). One describes 
the decay rate into an inclusive final state / in terms of the imaginary part of a 
forward scattering operator (the so-called transition operator) evaluated to second 



order in the weak interactions IfflL 72 



lmf(Q^f^Q) = Im / d 4 x %T (/V(x)£^(0)) (7.1) 

where T denotes the time ordered product and Cw is the relevant weak Lagrangian 
at the normalization point higher or about tuq. The space-time separation x in 
Eq. ( |7. 1| ) is fixed by the inverse energy release. If the latter is sufficiently large in 
the decay, one can express the non-local operator product in Eq. (|7. 1| ) as an infinite 
sum of local operators Oj of increasing dimensions. The width for Hq — > f is then 
obtained by averaging ImT over the heavy-flavor hadron Hq, 

(H Ql I m f(Q ^ f ^ Q)l H Q ) 

2M Hq v q JJ 
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r^W/ l 2 -m 5 -(/)/. a ( H Q I °i\ H Q)» (79) 
^ f \ v ckm\ rrig^Ci [fi) — — {(.2) 



2M Hq 

with Vckm denoting the appropriate combination of the CKM parameters. A few 
comments are in order to elucidate the content of Eq. ([7.2|) . 

(i) The parameter \i in Eq. (|7.2| ) is the normalization point, indicating that we 

explicitly evolved from mg down to fi. The effects of momenta below fi are lumped 

into the matrix elements of the operators Oi. 

( f) 

(ii) The coefficients c\ (//) are dimensionful, they contain powers of I/tuq that 
go up with the dimension of the operator Oj. Sometimes it is convenient to introduce 

dimensionless coefficients c| (ju) = mg"" 3 c- (ji) where di denotes the dimension of 

( f) 

the operator Oj. The dimensionless coefficients c\ depend on the ratio of final- 
to initial-state quark masses. Using the normalization introduced in Eq. ( [7.2| ), one 
obtains on dimensional grounds 

~(/)/ \ 1 ITT inlrr \ <n ( A QCD Cts^' 

with di denoting the dimension of operator Oi. The contribution from the lowest- 
dimensional operator obviously dominates in the limit mg — > oo. 

(iii) It seems natural then that the expansion of total rates can be given in powers 
of 1/rriQ. The master formula ( |7.2| ) holds for a host of different integrated heavy- 
flavor decays: semileptonic, nonleptonic and radiative transitions, CKM-favored or 
suppressed, etc. For semileptonic and nonleptonic decays, treated through order 
l/rriq, it takes the following form: 

G 2 F m%, 



M) ( J H q\QQ\ H q)m , C (f) (a]rr -2 (H Q \Q^GQ\H Q ) M 
2M Hq 5 W Q 2M Hq 



+E^MV <gQlwri 2 g g Q)|g ° >M + 0(17,4 



(7.3) 



We pause here to make a few explanatory remarks on this particular expression. 
First, the main statement of OPE is that there is no correction of order 1/mg 
||73|| . This is particularly noteworthy because the hadron masses, which control the 
phase space, do contain such a correction: Mh q = tuq (l + A/rriQ + 0(l/mg)); the 

parameter A, different for different hadrons, does not enter the width! The reason 
for the absence of the 1/mg correction in the total widths is two- fold: the corrections 
to the expectation value of the leading QCD operator QQ^ is only ~ /i 2 /mg, and 
there is no independent QCD operator of dimension 4 for forward matrix elements. 
Since the coefficients functions are purely short- distance, infrared effects neither can 
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penetrate into them. The absence of the dimension-4 operator in HQET was noted 
in 



72, 74 



A physically more illuminating way to realize the absence of corrections of order 
1/rriQ is to realize that the bound-state effects in the initial state (mass shifts, etc.) 
do generate corrections of order 1/mq to the total width - as does hadronization 
in the final state. Yet local color symmetry demands that they cancel against each 
other, as can explicitly be demonstrated in simple models. It is worth realizing 
that this is a peculiar feature of QCD interactions - other dynamical realizations of 
strong confining forces would, generally, destroy the exact cancelation. 

Second, the leading nonperturbative corrections are ~ (9(l/mg), i.e. small in 
the total decay rates for beauty hadrons: (/i/m b ) 2 ~ few % if // ^ 1 GeV. The first 
calculation of the leading nonperturbative corrections in the decays of heavy flavors 



was done in |73|, [75|, [76], p5| . 

Third, the four-quark operators (QTiq)(qTiQ) depend explicitly on the light- 
quark flavors denoted by q. They, therefore, generate differences in the weak tran- 
sition rates for the different hadrons of a given heavy flavor.^ Their effects were 



calculated already in mid-eighties W7 



Fourth, the short-distance coefficients (/j) in practice are calculated in pertur- 
bation theory. However it is quite conceivable that certain nonperturbative effects 
arise also in the short- distance regime. They are believed to be rather small in 
beauty decays f78j. 



Fifth, a new matrix element appearing in OPE, not discussed so far, is the scalar 
heavy quark density. Its expansion originally established in Ref. [nj takes the form 



(H Q \Q (it 2 + i<rG) Q\H Q ) 
(Hq\QQ\H q ) = (Hq\Q^ Q\Hq) + V 2 J 2 I +0{l/m%) . (7.4) 

Since (Hq\Qj Q\Hq) = 2Mh q , the spectator ansatz indeed emerges as the asymp- 
totic scenario universal for all types of hadrons, and holds up to 1/rriQ corrections. In 
addition to Qif 2 Q, the second dimension-five operator is the chromomagnetic oper- 
ator QioGQ. Since QD 2 Q is not a Lorentz scalar, it does not appear independently 
in Eq. Q. 

Equation ( [7.4| ) is readily obtained in the heavy quark expansion if one uses proper 
non-relativistic heavy-quark spinors incorporating the Foldy-Wouthuysen transfor- 
mation. In the context of HQET this was suggested in Ref. [jnj . When the Foldy- 
Wouthuysen transformation is ignored, the straightforward evaluation of the scalar 
density in HQET leads to incorrect l/m^ terms. Recovering the additional terms 
required a certain revision of the standard HQET strategy |79|, j80[| . At present, the 
full-QCD derivation and that based on HQET are in perfect agreement. We thus 

8 Expanding {Hq\Qi(jGQ\Hq) / tuq also yields contributions of order 1/mg; those are, however, 
practically insensitive to the light quark flavors. 
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a b 

Figure 4: a) Spectator effect leading to 1/rrib correction in the decay width b — > 
s + 7. b) A different cut of the same diagram leads to an electromagnetic correction 
in the hadronic decay width. Terms l/m& and cancel out in the sum of two 

decay widths. The solid dot in the vertex denotes the penguin-induced b — > sqq 
interaction. 



have, for example, for the pseudoscalar mesons, 

1 2 3 M 2 ikf 2 

2M Pq 2m 2 Q 8 m 2 Q 

The reason for the kinetic operator term to appear is quite transparent. The first 
two quantities on the right-hand side of the equation represent the mean value of 
the factor yl — v 2 reflecting the time dilation which slows down the decay of the 



quark Q moving inside Hn [pq, 81, 82 



Finally, Eqs. ( |7.3| )-( |7.4|) show that the two dimension-five operators do produce 
differences in B versus Ab/^b versus O5 decays of order 1/toq. To a small extent they 
can also differentiate B and B s via the SU(3) breaking in their expectation values. 
Differences in the transition rates inside the meson family are generated at order 
I/tuq by dimension- six four-quark operators. They are usually estimated in the 
vacuum saturation approximation which - although cannot be exact - represents a 
reasonable starting approximation. There is an intriguing way to check factorization 



experimentally ||81|| : similar four-fermion operators enter semileptonic b — > u transi- 
tion rates. Moreover, in the heavy quark limit the four-fermion operators populate 
mainly the transitions into the hadronic states with low energy, and thus show up, 
for example, in the end-point domain of the lepton spectrum where their relative 
effect is enhanced. Considering the difference of the decay characteristics of the 
charged and neutral 5's in this domain, one can measure these matrix elements and 
even feel their scale dependence. Further details regarding the "flavor- dependent" 



preasymptotic effects can be found in |83 



It is important to keep in mind that the OPE approach discussed in this sec- 
tion implies that all decay channels induced by a given term in the short- distance 
Lagrangian in Eq. ( [7.1[ ) are included. It is not enough to consider the states that 
appear at the free-quark level. The final-state interactions can annihilate, for exam- 
ple, the cc quark pair into light hadrons; electromagnetic interaction, if considered, 
can do the same. Disregarding the channels that can emerge due to such final state 
interactions can violate the general theorems: even I/thq terms can appear in such 
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incomplete "inclusive" widths. For example, the correction to the 6^5 + 7 width 
does have nonperturbative corrections scaling like 1/mj, due to the effect of weak 
annihilation (in mesons, or weak scattering in baryons) of the light quarks with 
emission of a hard photon in the penguin-induced weak decay b — > sqq, see Fig. 4a. 
This effect would only cancel against the (virtual) electromagnetic correction to the 
hadronic penguin- induced width, Fig. 4b. [|4], jSlfl . 



8 Applications of Heavy Quark Expansion 

If in the first part of the review we were more concerned with the theoretical founda- 
tions, here we pass to applied problems. The number of applications is quite large, 
and continues to grow. We will limit ourselves to those applications where the heavy 
quark expansion is exploited as a tool for extracting such fundamental parameter as 
\Vcb\, the element of the CKM matrix. Below we dwell on two methods often briefly 
called inclusive and exclusive approaches. 



8.1 \V c b\ from the total semileptonic B width 

The idea of extracting |V^,| from the inclusive semileptonic transition rate b —>■ civ 
is quite evident. If m& — > 00 the parton formula becomes exact. Comparing the 
decay rate in this limit to the theoretical expression one could determine |V^,|. In 
the real world, where is finite, analysis of various pre-asymptotic corrections, 
which became possible after the advent of the heavy quark expansion, is needed. 

Dimensional arguments tell us that the total decay width scales as the fifth power 
of mass, but not which mass it is. A priori possible choices would be the heavy flavor 
quark or hadron masses which agree to leading order, {Mh q -mg)/nig ~ (9(A/mg). 
However, in the case of beauty Mj, m 1.5m|. This difference, which would represent 
a power-suppressed nonperturbative correction, has to be brought under theoretical 
control if one wants to extract reliable numerical values for the CKM parameters 
from the measurement of the total semileptonic B width. In particular, the theo- 
retical uncertainties intrinsic to the pole masses must be exterminated. 

The heavy quark expansion yields an unambiguous answer. Not surprisingly, 



everything is expressed in terms of the quark masses |73| , |76| . [35 1 





where ellipses stand for higher order perturbative and/or power corrections; z , 
are known parton phase space factors depending on rv?Jm\ (see, e.g. |85|). The l/mf 
power corrections to r sl are rather small, about —5%, and thus the direct impact 
of the higher-order power corrections is negligible. One obviously has to be careful 



in dealing with the quark mass in Eq. ( |8.1|) since it enters in a high power. Our 
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strategy is to first present a complete theoretical formula which explicitly exhibits 
theoretical uncertainties. Then we inspect the extracted value for \V cb \. Finally, we 
detail the criticism that sometimes is raised against such an analysis. 

The first-order perturbative correction term iw? c jm^) is essentially the same 
as for fi — > evv known since the mid-fifties PB] . Recently the QCD calculation was 



improved by computing the second-order [87| and all-order |85j corrections arising 
due to the running of a s in the first-order diagrams (i.e. the BLM part of a k s 
corrections). These are presumably dominant terms. They have a small numerical 
impact on the width, though. For completeness, the second-order BLM terms are 
included below. 

Evaluating Eq. ( |3.1|) we find 

IVW = 0.0419 filial * {™*2L\ * ■ U - 0.012-<"* " 5 ^ 



\ 0.105 J V t b J \ ' 0.1 GeV 2 J 

f 1 - 0.01^4) (l + 0.006 a?5(1GeV) - a336 ) (l + 0.007^^) . (8.2) 
V 50MeV/ I 0.02 J \ 0.1 GeV 3 J V ; 

Here we have relied on the mass m&(/z), see Eqs. ( p.!4|) - (|3.15|) , and on the value 
of the low-energy coupling a^ s (lGeV) ~ ct^(2.3GeV) = 0.336 obtained in the 
dedicated analysis of pl| . The charmed quark mass m c (fi) is then obtained using 



the relation ( |3.16| ) which introduces also the parameter p. Furthermore, p 3 reflects 
the dependence on the 1/tuq terms. Electroweak corrections are still neglected. 
It is worth noting that all expressions for the perturbative coefficients depend on 
the concrete definition of mass and relevant operators. This dependence propagates 
into numerical expressions, which, however, give the same final result for observables 
for commensurate input. The numbers shown above correspond to the definitions 
discussed at length in the preceding sections. 

The proper evaluation of the theoretical uncertainties is a real challenge which 
we now discuss. The last four terms in parentheses in Eq. (|8.2|) exhibit some of the 
uncertainties; others will be considered shortly. 

(i) As was mentioned, the perturbative corrections are known exactly to first 
order in a s . Their impact on \V c b\ ~ (r^ 1 ^ J is only about 5%. The (presumably 



theor 

dominant) BLM part of higher-order corrections is now known to all orders; already 
in the second order the impact is less than 1%. The only remaining uncertainty, thus, 
is the genuine, non-BLM second-order corrections. They are not known completely, 
but there are good reasons to believe that they are indeed small. 

There exist some enhanced higher-order corrections not related to the running of 
a s , that are specific to the inclusive widths p9| . They appear if one uses the masses 
normalized at a high scale fi > m^. They have been accounted for in the analyses 
41, ESfl, but were not included in Refs. IBB, BHl. 



In the SV limit, when 

£ = (m b - m c )/{m b + m c ) < 1 
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the perturbative corrections to the widths can be related to the zero-recoil 
renormalization of the axial current, which was recently calculated by Czarnecki 
to two loops. In this way one obtains for the second-order non-BLM coefficient 



(8.3) 

where are defined in Eq. ( |3.9|) below. If ^ = it is clear that a°(C) is rather 
small numerically. At the actual value of £ ~ 0.5 the coefficient a® is different, of 
course, but hardly can be large if the appropriate value of p is used. For example, 
the second-order non-BLM term in dr s i/dg 2 at q 2 = g^ ax = (m,& — m c ) 2 constitutes 
approximately 0.1(a s /7r) 2 for the actual quark masses, according to Ref. The 
same correction to dr s i/dg 2 at q 2 = was recently evaluated ~ 1.25(ct s /7r) 2 in 



Ref. pi| , thus demonstrating that the non-BLM effects are moderate. Altogether, 
we assess a ±2% uncertainty in T s \ due to higher-order perturbative corrections, 
as a conservative estimate. There is an additional uncertainty associated with the 
experimental error bars in a s , but it is minor. 

(ii) The power corrections in r s i scaling like Aq CD / (m^— m c ) 3 , are at the 1% level, 
i.e. tiny. They emerge from several sources. First, there are l/m& corrections in the 
expectation values of the kinetic and chromomagnetic operators over the actual B 
meson state. They are small and, in any case, covered by the existing uncertainties 
in the values of the condensates. The "genuine" corrections are given by the local 



Darwin term p 3 D and the "spin-orbital" expectation value p 3 LS [[H]] appearing in 



Eq. (|2.11|) . The latter, in turn, appears also as a (local) 1/mj, correction in the 



expectation value of (B\b^aGb\B) which enters the 1/m 2 term in the widths. Thus, 
this effect is also negligible. It is worth noting that p 3 LS is expected to be suppressed 



in the ground-state mesons [10 



The largest direct effect is due to the Darwin term; relying on the approximation 
m 2 jm\ <C 1 and using the calculations of Refs. [^, ^(J one arrives at the estimate f\ 



r„(6 - c) oc L - O.Ol^^U = „ (l - 0.02^1^) . (8.4) 



Finally, if m c is related to m& via the mass formula fl3.16|) - so far, the most 
accurate method - one has an indirect dependence on p 3 D and p 3 through — m c . 
Even though the value of p 3 D is reasonably well estimated, the overall dependence 
on it practically cancels out in T s \ in this method. As a result, the extracted value 
of \V c b\ is sensitive only to non-local correlators p 3 n and p| as far as the 1/toq terms 
are concerned. Both p 3 n and p| are positive, see [HJ. The corresponding shift in 
T sl is thus, probably, negative but below 2%. 

9 There is a minor disagreement in the relevant expressions in the two works. Numerically it is 
negligible. 
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(iii) The calculation of the width could be affected by violations of duality, which 
are conceptually related to the asymptotic nature of the OPE expansions. This is the 
least understood ingredient of the theoretical analysis. Some model considerations 



are given in Ref. |[7ql . With the energy release in the semileptonic widths > 3.5 GeV, 
the duality violation is expected to be negligible, below 1% level, and not to exceed 
the effect of 1/ mn corrections. 

Assembling all pieces together and assuming, conservatively, that at present 
\Srrib\ < 50MeV and \p\ < 0.1 GeV 3 , we arrive at the model-independent estimate 

1^1 = 0.04,(25^)*^% 



2\ 



( (al - 0.5 GeV : , , 

I 1 " - 012 oTGeV^J ' ( X ± °- 015 P-t ± 0.01 mb ± 0.012) , (8.5) 

where the last error reflects nig 3 and higher power corrections, including possible 
deviations from duality. 

8.1.1 F(B -> XJv) 

Similar to the treatment of T(B — > X c lv) it is straightforward to relate the value of 



\V ub \ to the total semileptonic width Y{B — > X u lv) [28 



W4 = 0.00465 ( mB l; )2 XM y (^f • (1 ± 0.025 pcrt ± 0.03„„) . (8.6) 

The dependence on is practically absent here. 

Accurate measurement of the inclusive b — > uiv width is difficult and for a 
long time seemed unfeasible. However, recently ALEPH announced the first direct 
measurement p3fl : 

BR(fi -> XJv) = 0.0016 ± 0.0004 . 

It is not clear how reliable and model-independent are the quoted error bars in 
this complicated analysis. It certainly will be clarified in the future. Taking these 
numbers at their face value, one would arrive at the model-independent result 

\V ub \/\V cb \ = 0.098 ±0.013 . (8.7) 



The theoretical uncertainty in translating T(B — >• X u tv) into \V U \\ p8f is a few times 
smaller than the above experimental one. 

8.1.2 Caveat 

Deviation from duality is the most vulnerable element of the theoretical prediction. 
Reliable information on this aspect of QCD is scarce. The only model specifically 
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designed to address the issue is presented in Ref. [0. If it can be relied upon, at 
least for orientation, we can safely neglect deviations from duality in the semileptonic 
decays of B mesons, but not in those of D mesons. And indeed, a parallel analysis of 



Y(D — > X s lv) p2| |94|. p2| , along the lines discussed above, shows that the expression 
for T(D — ► Xgiv), similar to Eq. ( |8.1| ), falls short of the experimental number 
roughly by a factor of two, provided that the modern values of the fundamental 
parameters (m c , /x^ and so on) are used. 

8.1.3 Why has T S \(B) been sometimes discarded in the quest for |V^,|? 

The status of the radiative corrections to the widths had been the subject of some 
controversy recently, but finally the issue seems to be settled. The reason for appar- 
ently contradicting opinions was related, once again, to the problem of the heavy 
quark mass in the context of the heavy quark expansion. We will briefly sketch the 
problem and explain the answers. 

The inclusive width can be expressed in terms of any well-defined mass parame- 
ter. Its perturbative component through a certain order can conveniently be given in 
terms of the pole masses. Yet they are not defined with the necessary accuracy so as 
to include power-suppressed corrections. As discussed before, the pole mass contains 
long-distance contributions of order Aqcd- With the width being a short- distance 
quantity we can then infer that the indefiniteness ~ Aqcd i n the quark pole mass 
has its counterpart in an irreducible uncertainty ~ Aqcd/^q in the perturbative 
corrections if those are evaluated with pole masses as input. This conclusion has 
been proven in the most general terms by virtue of OPE where the emergence of the 
1/tbq infrared renormalon was noted ||22|| . It was also substantiated by the pertur- 



bative calculation. A concrete analysis in the framework of the BLM resummation 



was performed later in |5!| |85| . Consider the perturbative expansion of r s i in terms 
of the pole masses: 

GplTlf, pole,,, 1 2 / m cpole\ (si) (si) 

1 - ! — \Vcb\ ■ Z , " 4 > C; 



1927T 3 



, m b pole, 



3 5 °3 




One finds that ~ —1.7 is followed by a huge coefficient |87], ~ — 10 to 

—20 for b — > c (the exact value depends on the choice of the expansion parameter 
a s ) or even c4 — — 30 for b — > u. If this were the end of the story, this fact 
would mean that the a 2 s corrections by themselves reduce the width significantly, 
namely, by about 10%, suggesting, at first sight, that the unknown third (and higher) 
order corrections are sizeable. In that case one could argue that the theoretical 
uncertainties in extracting \V c b\ from T s i(B) are considerably larger than stated 
above. It would be natural to be concerned [B7[] whether such an extraction can be 
trusted at all! 

The key point that eliminates the concern, is as follows. When one calculates 

(si) 

the radiative correction factor c 3 through order k, one has to evaluate the pole 
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mass likewise to that order - since the pole mass is only a formal perturbative 
construction. Then one finds that receives sizeable increase as well; m po i e 

breaths as k changes. The shifts are correlated so that the width remains almost 
the same numerically. 

In doing so one must be careful always to use the very same perturbative ap- 
proximation both in masses m^ le and in the coefficient C3. Speaking theoretically, 
the perturbative approximations to the width computed in such a strictly correlated 
way contain no uncertainty ~ 0(Aqcd/"^q)- Uncertainties ~ (9(AQ CD /mg) still 
remain due to spurious terms ~ {Sm po i c / mg) 2 having no OPE interpretation. On 
the practical side, the latter problem emerges at the one-percent level only - unless 
one attempts to continue to too high orders. 

If one, however, exploits the short-distance masses, say the masses normalized 
at the scale ~ 1 GeV, all these problems are completely avoided and one does not 
have to rely on what seems to be a numerical miracle. Indeed, one can use a fixed 
number for m^ fi which, unlike m po i e , does not breath. Moreover, one finds then that 

Tsl = ^iw^ IKb|2 ' Zo (Sf^y) ' ( x + aM K/7r) + aM K/7r)2 + -) (8 - 9) 



with ai(jj) —1, a,2 (/i) ~ 1, etc. [p8| . The exact values of a/c(/i) depend on the 
scheme used. The same is true for mf, iC (/z). The scheme dependence affects the 
overall result for T s i only at a percent level. Such variations are unavoidable since 
the perturbative series are truncated. 

To reiterate, the suspicion of a large and uncontrollable impact of the perturba- 
tive corrections on the absolute values of the semileptonic widths, quite popular in 
the literature in the past, was mainly due to theoretical subtleties with the pole mass. 
In particular, it was tacitly assumed that the quark pole mass has an unambiguous 
value, exact up to a hundred MeV, or so. Then it was observed that: 

(a) It is difficult to accurately extract m^ ole from experiment. In any particular 
calculation one can identify effects left out, which can change its value by ~ 200 MeV. 
This uncertainty leads to a theoretical error 5i ~ 10% in T s i(B). 

(b) When routinely calculating T s i(i?) in terms of the pole masses, there are 
significant higher order corrections 5u ~ 10%. 

Thus the conclusion was made: T s i(i?) cannot be calculated with accuracy better 
than ~ 20%, and, correspondingly, at best ~ 10%. 

Both observations (a) and (b) above are correct, beyond any doubt. If one more 
step is taken, however, the conclusion is invalidated - one should take into account 
the fact that the origin of these two uncertainties is actually the same, and therefore 
they practically cancel each other. 

Let us mention that the problems associated with the the pole mass sometimes 
surface in a superficially different form, leading (in the past) to the same feeling, 
that 5\ Vcb\/\ V c b\ is as large as ~ 10%. Quite often, inflated error bars in m h are 
quoted, see e.g. [p6| , where the uncertainty in m;, is set 200 to 300 MeV. One may 
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notice that this variation actually reflects the change in the values of mj, emerging in 
different perturbative definitions of mj,. Such a procedure is like comparing a quan- 
tity renormalized in the MS and MS schemes at the same scale without including 
the non-trivial translation between those two definitions. The large perturbative 
uncertainty ~ 20% in the value of T s \ in was deduced not from the analysis 
of the absolute prediction for the semileptonic widths. Instead, it was based on a 
comparison of the perturbatively-improved calculation of T s \ to an ad-hoc quantity 
dubbed r tree ; the latter was supposed to represent a numerical estimate of T sl be- 
fore any radiative corrections are included. However, r tree involves the fifth power 



of the mass which was presumed to be m po i e in | 96fl . Then everything depends on 
what is the value of m po i e in the tree-level calculations. The large uncertainty stated 
in can be traced back to the assumption that this tree-level mass is given by 
the actual all-order sum of the perturbative series emerging when one attempts to 
calculate the would-be pole mass of the b quark in terms of fhb{rrib)- This procedure 
relies on summation of the series similar to Eq. ( |3.7|) , which is particularly unnatu- 
ral in the context of the tree-level evaluation. A finite value for this tree-level mass, 



~ 5.05 GeV, used in ||96[| , was due to a trick with a non-summable series, a principal- 
value prescription for the Borel integral, and so was the ratio r/r tree = 0.77 ± 0.05. 
Values other than 0.77 could have been obtained as well, invoking other prescrip- 
tions. In any case, it is clear that such definition of r tree is not appropriate for 
evaluating the impact of the perturbative corrections. 

The epilogue of the long story with the determination of \V c b\ from the inclusive 
width can be summarized as follows. After all refinements and thorough analysis of 
the corrections, we are basically back to square one. The expression for \Vd\ given 
in Ref. |44| has changed by less than one percent (if one inputs the same values of 
the experimental parameters as in WW). What became clearer in the last two years 
is the acceptable range of the underlying key QCD parameters entering the heavy 
quark expansion, and the size of theoretical corrections. 



8.2 \V cb \ from T{B D*tv) at Zero Recoil 

Introduction of the universal Isgur-Wise function |J was a crucial step in the evo- 
lution of heavy quark theory. One aspect is of a special significance in applications: 
the fact that it is normalized to unity at zero recoil @, [|, [|. As was emphasized 
later (see e.g. P7| ) one can exploit this feature for determinations of \V c b\ from the 
exclusive B — > D* semileptonic transition extrapolated to zero recoil. To this end 
one measures the differential rate, extrapolates to the point of zero recoil and gets 
the quantity |V^,Fd*(0)|, where Fp* is the axial B — > IvD* formfactor. To extract 
1 14b | it is necessary to know |Fd*(0)|, which, although close to unity, still deviates 
from unity due to various perturbative and nonperturbative corrections. In the real 
world 

F D .(0) = 1 + (—) + 5 1/m 2 + S 1/m s + ... = 



46 



l + o(—)+o(—=) +0 (-?—) +0(^1+... (8.10) 

V 7r / \ m c/ \m c m b / \ m 6/ 

The absence of I/tuq corrections in |-Fd*(0)| was noted in passing in Ref. 0. 
This fact was cast in the form of a theorem by Luke This is nothing but 

the heavy-quark analog of the Ademollo-Gatto theorem for the SU(3)& breaking 
effects which is routinely exploited in determinations of \V US \ from K — > niv and 
semileptonic hyperon decays. 

The task of precision determination of \V c b\ from the exclusive transition requires 



a detailed dynamical analysis of various pre-asymptotic corrections in Eq. ( |B.1(J| ). 
The perturbative part, albeit technically complicated, is at least conceptually trans- 
parent. The theory of 1/wt.q corrections is more challenging. 

The urgent need in evaluation of the 1/mn corrections in Eq. ( |S.10| ) for prac- 



tical purposes was realized quite early []IE |. In these days the theory of the power 
corrections in heavy quarks was immature, our knowledge was scarce, so that it 
was hard to decide even the sign of 5i/ m 2 - The general impression was that this 



deviation is quite small, ^ 2% 0]. At the present stage we learned much more 
about the nonperturbative and perturbative corrections. We discuss them in 

turn, and then proceed to numerical analysis and survey of the literature. At first, 
a brief excursion into the theory is undertaken to sketch the main ingredients of the 
direction where most of the advances have been achieved - the sum rules for heavy 



flavor transitions |44], [[(J. Although related to those presented in Sect. 4, these sum 
rules are not identical (with a few exceptions, see Sect. 8.2.2) since we account for 
I/tuq effects. They express the moments of the inclusive probabilities in terms of 
the quark masses and the expectation values of certain local operators. Since the 
transition probabilities are non-negative, the sum rules lead to constraints on the 
exclusive form factors. 

8.2.1 Heavy flavor sum rules: generalities 

The sum rules are derived in QCD using the standard methods of the short-distance 
expansion [|K|. One starts the analysis from the forward transition amplitude 



T^ 12 \q ;q) = ^jL J d 3 x dx e i( ^- iqoXo (B\iT(cT^b(x), bT {2) c(0))\B) (8.11) 

where I^ 1 ' 2 ) are some spin structures or, more generally, local operators. The tran- 
sition amplitude contains a lot of information about the decay probabilities. As 
usual in QCD, one cannot calculate it completely in the physical domain of q. The 
amplitude ( |3.11|) has several cuts corresponding to different physical processes [ff2 



The discontinuity at the physical cut q < Mb — y M|, + q 2 describes the inclusive 
decay probabilities at a given energy released into the final hadronic system. The 
cut continues further than the domain accessible in the actual decays, see Fig. 5. 
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Figure 5: Cuts of the transition amplitude in the complex go plane. The physical 
cut for the weak decay starts at g = Mb — (M|, + g 2 ) 1 / 2 and continues towards 
g = — oo. Other physical processes generate cuts starting near q = ±(Mj> c + q 2 ) 1 ^ 2 
(one pair); another pair of cuts originates at close values of go- An additional channel 
opens at g > 2rrib + m c . 



In QCD we calculate the amplitude Q8.11J) away from all its cuts. Essentially, it 
can be expanded in the inverse powers of the distance from the physical cut, e : 



e = M B - ^M 2 + q 2 - g = m b - yjm 2 + q 2 + 5(q 2 ) - g , (8.12) 

where 

5{q 2 ) = (M B -m b ) - {sjMl + <p- y/ml + q 2 ) ~ O ( A QCd£^) + O . 
The "off-shellness" e must be chosen in such a way that that |e| and |e-arg e| ^> Aqcd 



but, simultaneously, |e| <^ y m 2 + q 2 , m&. The second requirement allows us to 
"resolve" the contributions of the separate cuts. 

How does the deep Euclidean expansion of T^ 12 ) help constrain the amplitude 
in the physical domain of real e lying just on the cut? The dispersion relation 
immediately tells us that the coefficients in the expansion of the amplitude in powers 
of 1/e are given by the corresponding moments of the spectral density in T^ l2 \ 

T( If MT^Jl & , = _1 I ^ I ImT<121( , ?v * d( , , 

7r J e — e 7r , n e k+L J 

(8.13) 

Various subtleties going far beyond the scope of the present review, are discussed in 
Refs. and @. 

On the other hand, we can build the large-e expansion of the transition amplitude 
per se, treating it as the propagation of the virtual heavy quark submerged into a 
soft medium. The expansion takes the general form 

T (12) (»;9) = 1 l-Jd\dx ls e i ""»y 
(m c + i p- ^_ (;C ^ )2 _. ffG r«K0) W = 
48 



mj, - (m b - g ) + q — 2m b n + 2qn - n 2 - |crG 
(r«(m c + (E c + e - <5(<f)) 7 o + 97+ tf) x 

oo [2(E C + 6 - ggjjTTQ + 2g + IT* + | 
n=0 [-( e -5(g 2 ))(2E c + e -5(g2 ))] (n + D 



Here _E C = Jm 2 , + q 2 , and we use the short-hand notation 



(•■•) = J d*xdx e^-™*° (B\b(x)...b(0)\B) ■ 

In the first equation we used the full QCD fields and then passed to the low- 
energy fields according to Eqs. (|2.3|) - (|2.5|) . Picking up the corresponding term 1 / e k+1 
in the expansion over 1/e at e > Aqcd and evaluating the expectation value of the 
resulting local operators (e.g., &(0)7r M 6(x) = 8 A (x) &7r M &(0)), one gets the sum rules 
sought for. Taking k = yields the sum rule for the equal-time commutator of the 
currents cT^b and bT^c; k = 1 selects the commutator for the time derivative, etc. 
After this brief general overview we now outline a few concrete applications [TO . 



8.2.2 Sum rules for /i 2 and n G 

Considering the semileptonic transitions driven by the pseudoscalar weak current 
J5 — S d 3 x (ci^b) (x) (i.e., at zero recoil, q — 0) one obtains the sum rule for the 
structure function 

^ rw (5 \e)de = (l/2m c -l/2m b ) 2 L 2 M-ii 2 M) • (8-15) 
In Jo v / 

The normalization point /1 of fi 2 and \^ G is introduced through the cut-off of the 
integral on the left-hand side. Since the structure functions are non-negative, one 
arrives at the conclusion that a4(//) > A* 2 ? (/•*)) which is the field-theoretic analog of 
the difference of the sum rules ( |4.15|) and ( |4.14| ). 

The sum rule ( 4.14 ) for fi G can be also easily obtained in this way. For example, 



one can consider the antisymmetric part (with respect to i and j) of the correlator 
of the vector currents (T^ = ji, = jj). Another possibility is to turn directly to 
the sum rule for the correlator of the non-relativistic currents CKjb and banc in the 
leading order in 1/rriQ. The OPE guarantees that all such relations are equivalent. 

10 Corrections in the overall normalization of the current are ignored as irrelevant. The structure 
functions are defined, e.g. in |99|. 
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8.2.3 F D * at zero recoil 



The axial current 074756 produces Z?*, Dtt and higher excitations in semileptonic B 
decays at zero recoil. A straightforward derivation yields the following sum rule for 
this current, 



\F D . 
where 



+ 



2tt 



w A (e) de 



£>0 



Ufa) - A i 



/m 2 



ITT 



/m 2 



j4_ 

3m 2 . 



+ 



_L_ J_ 2 

x m 2 ml 3m c mb / 



(8.16) 



(8.17) 



and (D*(q = 0)|c7i75&|.£>) = \J2MbMb*Fb*^%- Furthermore, w A denotes the cor- 
responding structure function of the heavy hadron excited by the axial current. 
Its integral describes the contributions from excited charm states with mass = 
M D * + €i, up to 6j < \i: 



1 

2^ 



w A (e) 



de 



e>0 



£i<H 



Contributions from excitations with e higher than \i are dual to perturbative con- 
tributions and get lumped into the coefficient Ufa) of the unit operator, the first 
term on the right-hand side of Eq. (|8.16 ) . 

The role of \x is thus two-fold: in the left-hand side it acts as an ultraviolet 
cutoff in the effective low-energy theory, and by the same token determines the 
normalization point for the local operators; simultaneously, it defines the infrared 
cutoff in the Wilson coefficients. 

Equation flS.16|) immediately implies 

\F D *\ 2 = Ufa) - Af /m2 - £ \F, t \ 2 > Ufa) - Km 2 (8-19) 

due to the positivity of \Fi\ 2 . (For a moment we forget about the cubic corrections, 
A^ m 3 .) The perturbative coefficient Ufa) * s obtained considering the sum rules 
( |S.16| ) and (|8.15| ) in perturbation theory. To order a s []10 



Ufa) 



1+ 2— 

7T 



m b + m c m h 8 1 / /T u? 
In + - \ —^ H o + 



2fi 2 



mt 



mt 3m c m b 



(8.20) 



The most complicated part of the genuine second-order correction in Ufa) was 
calculated in Refs. [pop . It turned out to be small. The two-loop calculation of Ufa) 
was recently completed | |100| . The BLM-resummation of the one-loop result was 
carried out in |51| and placed the numerical value of Ufa) somewhere in between 
the tree-level and one-loop estimates, see Fig. 6. This analysis suggests that at 
fj, ~ 0.5 GeV, the value of fjAfa) = (Ufa)) 1 is 0.99 ± 0.01 with the uncertainty 
coming from the higher-order corrections and terms ~ (fi/m c ) 3 . 
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Figure 6: The value of rj A {fi) for Aqq D = 90 MeV, 130 MeV and 175 MeV. ^(/i) 
must be used in the Q CD-based calculations of the zero- recoil B — > D* formf actor 
when nonperturbative effects are addressed. 



The nonperturbative l/uin power corrections 6y m 2 to Fr>*, Eq. ( |8.10| ), are neg- 
ative. They consist of two terms: the first, — |A^ m2 , is known explicitly. The 
second, the sum over the excited states, also scales as 1/toq, it is negative, but is 
unknown otherwise. It cannot be calculated without additional dynamical assump- 
tions. — Sy m 2 exceeds ^/(6m^) ~ 0.035 in magnitude. The sum over the excited 
states replaces the wavefunction overlap deficit of the simple quantum-mechanical 
analysis of the transition probabilities. We will return to this, the most difficult part 
of the analysis, after briefly discussing 1/tuq terms in the heavy quark expansion. 

It is not difficult to calculate the I/tuq term in the sum rule. Its explicit form 
depends on the convention for the D = 5 operators in Eq. ( |S.17| ). Assuming that 



\x\ and [Aq there are understood as the expectation values of the operators b{iD ) 2 b 
and b ^aG b over the actual B mesons rather than the asymptotic states, we get 



n 3 1 „3 



Af = pD Ifag + __ ( _L _ __ ) {p l + p\ ) . (8.21) 

1/m 4m 3 c 6m c m b \m c 2m b J v D Lbl K ' 

If the asymptotic values of the D = 5 matrix elements are used, the additional 1 / mf 
parts (see come from the expectation values in Eq. ( |8.17 ). The l/rnq correction 



to the inclusive sum of the decay probabilities turns out to be moderate, and we 
will not dwell on this issue further. 

The excited states in the sum rule are higher excitations of D* and non-resonant 
Dtt,... states. Their contribution in the right-hand side of Eq. ( 8.19Q is negative. 



What can be said about its absolute value? Unfortunately, no model-independent 
answer to this question exists at present. The best we can do today is to assume 
that the sum over the excited states (cut-off at //) is a fraction x of the local term 
given by /4 and fy, 



1 



w A {e)de = £ = X A? /m2 , (8.22) 



27r Joo a<n 
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where on general grounds x ~ 1- Q Some idea of the value of x can be obtained from 
the Dir contribution which is readily calculable by virtue of the soft-pion technique. 
The result exhibits the proper scaling behavior in all relevant parameters Its 



absolute value depends on the pion coupling to the heavy mesons. The estimate 
of this constant was improved |101| . With this improvement, the Dir contribution 



to x fells i n t ne interval 0.15 to 0.35. It is natural to think that the resonant 
contribution is several times larger. To be conservative for numerical estimates we 
accept x = 0.5 ±0.5. Rather arbitrarily we limit x by unity on the upper side. Note 
that the larger the value of x the stronger the deviation of Pd*(0) from unity. In 
this way we arrive at 



D* 



Va{h) ~ (1 + x) 



6m?. 8 \ml ml ?>m c m\, 



5x/ m 3 . (8.23) 



Assembling all pieces together and assuming our standard input for /i 2 and the 
low-scale quark masses (m c = 1.25 GeV, m b = 4.75 GeV) we get 

Fd* ^ 0.91 - 0.013 ^"^^ ± 0.02 cxcit ± 0.01 pcrt ± 0.025 1/m 3 . (8.24) 

Estimates of the uncertainties in the 0{l/mq) corrections and the contributions 
from the higher excitations are not very firm and reflect a reasonably optimistic 
viewpoint of the convergence of the heavy quark expansion. Altogether we get 



D* 



0.91 ±0.06, (8.25) 



where the optimistic uncertainty ±0.1 GeV 2 is ascribed to /x 2 . 

The uncertainty can presumably be reduced in the future by one percentage point 
by accurately measuring /i 2 . An additional 0.01 can be removed via a dedicated 
experimental study of the transitions to the excited P-wave states. Using current 
theoretical technologies, it seems impossible to overcome the 5% barrier of the model- 
independent theoretical accuracy in this exclusive decay. 

In summary, the QCD-based analysis favors a significantly larger deviation of 
Pd*(0) from unity than those en vogue three years ago. The 1/m 2 . shift by —0.035 is 
definitely model-independent. Values of \Si/ m 2 1 close to the minimal possible scenario 
0.035 are achievable only if there are practically no zero-recoil transitions B — > 
Dtt, D**, ... up to the mass range 2.5 to 2.6 GeV, and, additionally, the overall yield 
of the lowest P-wave states in the sum rule for is totally suppressed. Needless 
to say, that this is hardly possible. In any case such a distinct pattern can also be 
checked experimentally. Although the present data are not conclusive yet, there is 
no evidence for such a feature. 



11 In perturbation theory \ = 1 in the first order, but this relation changes in higher orders. 
12 The central numerical value here is larger by 0.015 than that cited in p4| for the same jj^. 



0.01 is due to a shift in the perturbative factor rjA{n) (0.98 was used in 44 ). The remaining 0.005 
reflects a different definition of ^ at the level of perturbative corrections. 
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8.2.4 Quantum-mechanical interpretation 



The sum rules (pTTB) have a transparent interpretation in the framework of con- 



ventional quantum mechanics fllOfl . From the gluon point of view the semileptonic 



decay of the b quark is an instantaneous replacement of b by c quark. The over- 
all probability of the produced state to hadronize to some final state is exactly 
unity, which corresponds to the first term on the right-hand side of Eq. ( |8.16| ). The 
nonperturbative corrections in the sum rule appear since the normalization of the 
weak current 07^756 is not unity and depends, in particular, on the external gluon 
field. Expressing the current in terms of the non-relativistic fields used in quantum 
mechanics one has, for example: 

c 7fc75& ^a k -\ + a ^ dl ^ 2 _ (°iD)vk(aiD) \ + (—\ (8 26) 

8m;? 8m? 4m c mft / \m 3 



The second term immediately yields the correction to the normalization, which is 
present on the right-hand side of the sum rule. It is curious to note that the first 
HQET analysis of 1/m 2 corrections (see, e.g., [16|, |65| ) does not include the first 



two terms in the brackets which contain the dominant effect ~ l/ m c- They appear, 
however, in the approach of the Mainz group . 



8.2.5 Digression in the literature 

In view of the practical importance of the value of Fjj* (0) a more detailed comment 
on the literature seems in order. We will dwell on a combined analysis ||102|| often 



cited in this context as the analysis with the best accuracy and reliability available 
today on the theoretical market. Our task is to warn the potential reader of hidden 
assumptions and inconsistencies. 

The basic idea of Ref. | |102| | was supplementing the sum rules of [33, |ItJ by a 



certain symmetry relation obtained previously ||16|| . In this way one gets rid of the 
sum over the excited states for the price of expressing the deviation of the formfactor 
Fp* from unity in terms of two other unknown zero-recoil formfactors for the vector 
current in pseudoscalar and vector mesons £py, and an unknown hadronic correlator 



\ 2 



<*> = -[^r-^r)(^r-^r) + t^~t (~U + 6/4 - *a ) ■ (8.27) 



1//m \2m c 2m;,/ \2m c 2m^ J Amcm^ V 3 

Here lp and ly are defined via the forward matrix elements 

1 / / 1 1 

:{D\cj b\B) = r) V 1 



2^M B M D N 1 ' 1 ' ' \ \2m c 2m b 
1 . nl/ . u Ilm/ .,. , f ( 1 1 



2 



53 



(£ p is, theoretically, observable in the decays B — > r tv). The same sum rules derived 



in 



44], |10| for F D ( t ) apply to Ipy as well. The parameter X G is a sum of three other 



unknown correlators introduced in [|16j.0 The strategy of Ref. ||102|| was to obtain 



constraints on the introduced unknown parameters following from the sum rules 
( ^.16| )-( p.l5| ), and, then, in turn, to find the allowed interval for 5 1 / m 2. 

Expressing the phenomenologically relevant formfactor Fp* in terms of two other 
unknown formfactors and one unknown correlator is clearly of no help unless a new 
dynamical input is provided. It was provided in a form of an ad hoc prescription for 
the value of Xq which was set to be very small. This step was crucial since without 
information on X G the relation ( |8.27| ) provides no constraint. 

Another crucial element of the analysis was exploiting the quark model to esti- 
mate l P and ly. The quantities 



1 



1 



2m c 2m;, 



"p,v 



were identified with the deviations from unity of the wavefunction overlap between 
the charm and beauty states, associated with the c, b quark mass splitting. The 
complete relation has, however, the form [|1(| 



[*d\*b) = 1 



1 



1 



2m c 2mi 



+ 



a4 



5ij 



1 - 



2m r 2mt 
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where ^ B (*) D (*) are the normalized wavefunctions for charm and beauty. The terms 
with /i^ and p? G were missed in Ref. ||102|| . These terms emerge due to the Foldy- 



Wouthuysen transformation. Thus, even accepting the quark model for guidance, 
we observe inconsistency in evaluating the dominant 1/toq part of the correction. 
This led to the surprising statement on the insensitivity of the estimates for 5 1 / ?n 2 to 
the value of kinetic operator. This insensitivity, together with discarding 1 / rriQ cor- 
rections, was the reason behind a smaller overall uncertainty in Fjj* than indicated 
in Eq. ( ggg ). 



Another potential problem of this work refers to the treatment of radiative 



a" 



corrections. The problem arises once one addresses the issue of higher-order sum- 
mation ignoring the basic aspect of OPE. Postulating that the perturbative factor 
tja{^) in Eq. ( ^.23[ ) and its vector counterpart rjv(fi) coincide with the sum of purely 
perturbative corrections involving all virtual momenta (including those below Aqcd) 
we are bound to run into troubles. Physically, this is double-counting since the do- 
main below fi was already included in the matrix elements. It is explicit in this 



approach: the sum rules of |44|, [IIJ are valid only if the perturbative coefficients are 

13 As a matter of fact, its meaning is seen from relation ( |8.16| ) in the quantum-mechanical inter- 
pretation at mb — fn c . It represents a sum of the transition probabilities in the spin- flip processes. 
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understood in the Wilsonian sense; they definitely do not hold in the interpretation 
used in | 102|| . Formally the resummed series itself tells us of an inconsistency: the 



resummed values of the perturbative factors r/y, ?7a cannot be defined at the level 
(9(AQ CD /m^); they develop unphysical imaginary parts and so on. This irreducible 
ambiguity in thus defined t/a was estimated to be 2.4% | |103| . It is impossible to get 



rid of it in the analysis [[L02 |, since the hadronic parameters £p, iy etc. are treated 



as fixed positive numbers. This fact was pointed out in ||103|| and admitted in ||102 . 

A practical manifestation of these theoretical problems was later noted in Ref . ||104|| . 
It was found that already in the second order the incorporation of the perturbation 
theory in the approach of [ |102|| suffers from large perturbative corrections ~ 3-5%. 



Combining these observations it is logical to conclude that the combined analysis 
|102| | loses information rather than adds it and worsens the accuracy. It seems 



unlikely that following this line of reasoning in the future we could narrow the 
possible interval for Fp*. Evaluating deviations from the symmetry limit for lp and 
iy is neither easier nor more difficult than evaluating Fd* itself. At the moment 
there seems to be no ground for more optimistic attitude. 

8.2.6 Analyticity and unitarity constraints on the formfactor 

Extracting \ V c b\ from the exclusive decays implies extrapolating data to zero recoil. 
Near zero recoil statistics in the decays B — > D*iv is very limited, and the result 
for the differential decay rate at this point is sensitive to the way one extrapolates 
the experimental data to q = 0. Most simply this is done through linear extrapola- 
tion. Noticeable curvature of the formfactor would change the experimental value 
for \V c bFr>*(0)\. Under the circumstances it is natural to try to get independent 
theoretical information on the q 2 behavior of the formfactor near zero recoil. 

Some time ago it was emphasized ||105|| that additional constraints on the q 2 



behavior follow from analytic properties of the B — > D* formfactor considered as a 
function of the momentum transfer q 2 , combined with certain unitarity bounds. In 
particular, in the dispersion integral 

= I f ds ^4 (8.30) 



7T J s — q 2 

the contribution of the physical s-channel domain q 2 > (Mb + M D (*)) 2 is bound 
since \F\ 2 describes the exclusive production of BD and cannot exceed the total be 
cross section. The integral QS.30] ) receives important contribution from the domain 
below the open be threshold, where F(q 2 ) has narrow pole-like singularities corre- 
sponding to a few lowest be bound states. Introducing a set of unknown residues and 
making plausible estimates of the positions of these bound states, on the one hand, 
and adding a small non-resonant subthreshold contribution, on the other hand, 



Refs. ||106|| suggested an ansatz for the formfactor in the whole decay domain. 
Since the residues are unknown and the momentum transfer in the actual decays 
B — > D^*'iv varies in a rather narrow range, the advantages of this parameterization 
over the standard polynomial fit are not clear at the moment. 
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Surprisingly, a much more stringent (compared to what follows from Ref. [|106j| ) 
relation between the slope and the curvature of the formfactor (i.e., F'(q 2 ) and 
F"(q 2 ) at zero recoil) was found in [ |107| |. This relation was hastily incorporated in 
some experimental analyses. Below we argue that one should be extremely cautious 
in relying on the relation [ |107| ] due to a hidden assumption which makes the results 
very vulnerable and unstable. 

The suggestion of Ref. ||107f| was two-fold: consider the scalar formfactor for the 
B — ► D transition, and discard the resonant subthreshold contribution. In the heavy 
quark limit the scalar formfactor is related to usual axial- vector one F(q 2 ). If one 
neglects the contribution coming from the resonances below the open be threshold 
(in this channel they are scalar be relatives of the Xb an d Xc states) and retains only 
the the non-resonant continuum similar to the one used in ||106|| , stringent bounds 
do indeed follow from the dispersion representation. The non-resonant continuum, 
however, can be safely disregarded numerically. It is 50 to 100 times smaller than 

To give 



what can be expected from the resonant subthreshold contribution 
an idea why it happens it is sufficient to recall that in the good old pole- dominance 
models, the subthreshold poles due to Xbc saturate the formfactor completely. Rein- 
stating the Xbc subthreshold resonances at appreciable level suppresses the predictive 
power of the dispersion approach, as is explained in detail e.g. in Refs. | |109| |; no 
advantage compared to a more careful analysis of |pL 06|| emerges. 

An additional problem one immediately encounters is the necessity of including 
l/m c corrections translating the results for the scalar current into those referring to 
the axial formfactor. Such corrections to the derivatives over the velocity transfer 
were evaluated | 110 |, and turned out to be very significant even for more stable 
inclusive transitions. For further comments see Ref. 
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The lesson we would like to draw is as follows: it is not advisable to impose the 
model-dependent relation ||107|| in experimental extrapolations to zero recoil, where 
model dependence is undesirable. 



9 Challenges in Nonleptonic Beauty Decays 

Concluding the review it is impossible to avoid mentioning some unsolved mysteries 
of the heavy quark theory. Today's question marks carry the seeds of tomorrow's 
advances. Basically there are two problems where our theoretical understanding is 
lagging behind. Both are related to non-leptonic decays. The most placid solution 
(which we would prefer, of course) is experiment evolving towards theory, to meet 
its expectations. If this does not happen, some major adjustments in our theoretical 
ideas seem to be inevitable. 
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9.1 Semileptonic branching ratio of the B meson and T{B — > 

ccsq) 

The theoretical attitude to this problem oscillates with time. Twenty years ago it 
was believed that the parton model would give a sufficiently accurate prediction, 
BR(6 — > c£v) ~ 15%. Then it was realized that the a s radiative corrections were 
important The present situation is given in Ref. ||112|| . The general tendency 



associated with the a s corrections is lowering the value of BR(6 — > civ) down to 
~ 11.5 to 13.5%. This is still noticeably higher than the experimental number ||113|| : 

BR(E -+Xiv) = 10.43 ± 0.24% . (8.31) 

Nonperturbative corrections lower the semileptonic branching ratio further; yet be- 
ing of order \jm\ they are numerically quite small, A nonpert BR s i(i?) ~ 0.5%, which 



seems to be insufficient to close the gap between the expectation and the data ||114 

The issue of the semileptonic branching ratio must be considered in conjunction 
with the charm yield n c , the number of charm states emerging from B decays. To 
measure n c one assigns charm multiplicity one to D, D s , A c and S c and two to 
charmonia. Zero is assigned to the charmless hadronic final state. It is obvious that 

n c ~ 1 + BR(B ccsq) - BR(B no charm) . (8.32) 



The experimental situation with n c is as follows. The CLEO group finds n c 
1.134 ± 0.043 whereas the ALEPH collaboration reports n c = 1.23 ± 0.07. While 
both experimental numbers [|115|1 are consistent with each other and consistent with 



a "canonical" value 1.15, the CLEO number clearly favors lower values (~ 1.15) 
while ALEPH does not rule out n c as large as 1.3. The statistical average is n c = 
1.16 ±0.04. 

Further reduction of the theoretical prediction for the branching ratio can be 
provided either by higher-order perturbative corrections (which are amenable to 
analysis, at least, in principle) or by largely uncontrollable duality violations (usual 
scape-goat). In the later case it is natural to suspect T(B — > ccsq) since the energy 
release in b — > ccs is not very large. 

With the advent of the heavy quark theory the question of the compatibility of 
the existing theoretical ideas with the data on BR(I? — > X c lv) and n c acquired a 



solid footing ||1 14|| . While n c m 1.15 came out naturally, the excess of BR(I? — > X c £v) 



was an obvious challenge. Shortly after, the attitude changed. A natural desire to 



have all problems peacefully settled prevailed, see e.g. the summary talk where 
the general conclusion leans towards the absence of any problem. A combination 
of two factors was crucial in this respect. Large values of a s fashionable two or 
three years ago (corresponding to a s (Mz) = 0.125 or even higher) enhance the non- 
leptonic width and, hence, suppress BR(£> — > X c £v) down to 11 or even 10.5%. 
Simultaneously n c jumps up to 0.125, but since the ALEPH data were newer, their 
significance was overemphasized, and it was tempting to close one's eyes on the 
CLEO data. 
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The present perception of the issue seems more balanced; it acknowledges the 
existence of the problem, see e.g. 



116, 117 



The only theoretical ingredient one 
needs to reveal the problem is the statement that the b — > ccs channel is to blame for 
the discrepancy in the semileptonic branching ratio. It is conceivable that BR(5 — > 
ccsq) is actually larger than it is usually inferred from the explicit quark-gluon 
calculation, either due to higher order a s corrections or due to deviations from 
duality. Say, if BK(B — > ccsq) ~ 0.25, rather than 0.15, (which is in line with 
the fresh CLEO data) and the charmless modes are negligible, this would bring 
the predicted semileptonic branching ratio pretty close to the observed one, but 
n c becomes at least two standard deviations higher than the current average. A 
possible way out is to assume the charmless modes (say, b — > s + gluon) at the level 
of 10%. According to Eq. ( |8.32|) , this will bring n c down to the acceptable value. 
Reference ||116|| suggests that such high yields of b —>■ s + gluon are attainable in 

11711 new 



H suggests that such high yields of b 
QCD, which is extremely unlikely, to put it mildly. According to Ref. 
physics is supposed to be responsible. 



9.2 Lifetimes of Heavy-Flavor Hadrons 

As stated before, differences between meson and baryon decay widths arise already 
in order l/rug. The lifetimes of the various mesons get differentiated effectively first 
in order l/m%. A detailed review can be found in | 8~3fl ; here we will comment only 
briefly on the issue. 

Because the charm quark mass is not much larger than typical hadronic scales 
one can expect to make only semi-quantitative predictions on the charm lifetimes, 
in particular for the charm baryons. The agreement of the predictions with the 
data is surprisingly good. (It is quite possible, though, that future more precise 
measurements of the S c and fl c lifetimes might reveal serious deficiencies.) 

As far as the beauty lifetimes are concerned there is much less "plausible deniabil- 
ity" when predictions fail. Table [l] contains the world averages of published data Q 
together with the predictions. The latter were actually made before data (or data 
of comparable sensitivity) became available. 

Data and predictions on the meson lifetimes are completely and non-trivially 
consistent. Yet even so, a comment is in order for proper orientation. The numerical 
prediction is based on the assumption of factorization at a typical hadronic scale 
which is commonly taken as the one where a s (/^hadr) — 1- The related uncertainty 
has been emphasized in [|119|1 , however the recent QCD sum rule estimates ||120[| 



(though carried out in a simplified manner) did not support the conjecture of ||119 
about the significant impact of nonfactorizable contributions. While there is no 
justification for factorization at fi ~ m&, there exists ample circumstantial evidence 
in favor of approximate factorization at a typical hadronic scale - from the QCD 
sum rule calculations, to lattice evaluations, to l/N c arguments. More to the point, 



14 



It is unlikely that significant new data will appear before the turn of the millennium. 
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Table 1: QCD Predictions for Beauty Lifetimes 



Observable 


QCD Expectations (1/mb expansion) 


Ref. 


Data from B3 


r(B-)/r(B d ) 


l + 0.05(/V200MeV) 2 


H 


1.04 ±0.04 


f(B s )/r(B d ) 


1 ±0(0.01) 


||1 18 1 


0.97 ±0.05 


r{K b )/r{B d ) 


>0.9 


|118 


0.77 ±0.05 



the validity of factorization can be probed in semileptonic decays of B mesons in an 
independent way, as was pointed out in |8l| . 

The prediction on t(Aj ) ) versus r(B d ) seems to be in conflict with the data. 
However, the experimental situation has not been fully settled yet. The difference 
between (T(Ab)/r(£>d)) cxp . ~ 0.77 and t(Ai,) /T(B d )\ t \ lcov ~ 0.9 represents a large 
discrepancy. A failure of that proportion cannot be rectified unless one adopts a new 
paradigm in evaluating baryonic expectation values of the four-fermion operators. 
Two recent papers [|121| , |119|| have re-analyzed the relevant quark model calculations 



and found: 

r(A b )/T(B d ) = 1 - A (AbB) , A(a 6 _b) ~ 0.03 to 0.12 . (8.33) 

There are large theoretical uncertainties in A(\ b B) since the baryon lifetimes reflect 
the interplay of several contributions of different signs. Yet one cannot boost the 
size of A(a 6j b) much beyond the 10% level: to achieve the latter one had to go beyond 
a usual description of baryons when light quarks are "soft". A similar conclusion 
has been reached by the authors of Ref. | |122| | who analyzed the relevant baryonic 
matrix elements through QCD sum rules. 



10 Conclusions and Outlook 

Heavy quark theory is now a mature branch of QCD. Many practically important ap- 
plied problems that defied theoretical understanding for years, are now tractable. At 
the same time, all natural limitations of QCD take place also in the heavy quark the- 
ory. The infrared part of dynamics is parameterized rather than solved. Therefore, 
every new success based on the general properties of the quark-gluon interactions, 
is a precious asset. The most important stages of the success story are heavy quark 
symmetry itself, introduction of the universal Isgur-Wise function and combining 
Wilson's approach with the heavy quark expansion. The exact inequalities of the 
heavy quark theory is another link of the same chain. 

The most clear-cut recent manifestation of the power of the heavy quark theory is 
the framework it provides the determination of |V c b|. It is remarkable that the values 
of \V c b\ that emerged from exploiting two theoretically complementary approaches 
are very close. The progress was not for free: it became possible only due to essential 
refinements of the theoretical tools in the last several years, which prompted us, 
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in particular, that the zero-recoil B — > D* formfactor Fp* is probably close to 
0.9, significantly lower than previous expectations. The decrease in Fp* and more 
accurate experimental data which became available shortly after, reduced the gap 
between the exclusive and inclusive determinations of \V c b\. There is a hint in the 
experimental data that some discrepancy may still persist: the central value of \V c b\ 
from B — ► D* decay seems to be somewhat lower than that from T S \(B). Both 
theoretical values, however, depend to a certain extent on the precise magnitude of 
/i 2 , as is seen from Eqs. (|8.5| ) and ( |8.24| ). In the exclusive and inclusive formulae the 
dependence is rather similar in magnitude but opposite in sign. It is tempting to 
think that the actual value of /i 2 is somewhat larger than the "canonical" 0.5 GeV 2 . 
Increasing it by about 0.2 GeV 2 makes the two results much closer. We hasten to 
add, though, that the existing experimental error bars are such that any speculations 
on /x 2 are premature. Moreover, theoretical uncertainties in the exclusive formfactor 
also preclude us from the above adjustment of /i 2 . (Fd* can well be, say, 0.87 even at 
the canonical value /i 2 = 0.5 GeV 2 ). Future accurate measurements will, hopefully, 
allow one to directly measure - through comparison with T s i(B) - the exclusive 
formfactor with accuracy better than that achieved by today's theory. Thus, we 
will get a new source of information on intricacies of the strong dynamics in a so far 
rather poorly known regime. 

A large number of applications of the heavy quark theory are based on duality. 
Although this notion becomes exact at asymptotically high energies, at finite ener- 
gies (momentum transfers) certain deviations must be present. How fast duality sets 
in and how large are these deviations are important questions. These and similar 
questions are among most difficult, with virtually no or very little progress. De- 
terminations of | V c b | we discussed rely on the assumption that approximate duality 
between the actual hadronic amplitudes and the quark-gluon ones sets in already at 
the excitation energies ~ 0.7 to 1 GeV. While there are no experimental indications 
so far that this is not the case (at least, in the semileptonic physics), the proof is 
not known either. If that is not true, and duality starts only above 1 GeV, most 
probably one would have to abandon the idea of accurate determination of \V c b\ from 
the exclusive B — > transitions. The only option still open will be the inclusive 
semileptonic decays where the energy release is large, ~ 3.5 GeV. Of course, in such 
a pessimistic scenario (which, we believe, is unlikely) the theoretical precision in 
\V c b\ will hardly exceed 5%. 

What lies ahead? There are problems (e.g. duality violations) which we simply 
do not know how to attack. Solution of other problems seem to be possible in the 
near future. A practically important problem of this type is perturbation theory in 
the context of Wilson's approach, where the soft parts of all diagrams have to be 
removed from the a s series. This will lead us to more accurate estimates obtained 
in the heavy quark expansion. Another topical problem is constructing a reference 
model of the semileptonic decays where the transition amplitudes will be saturated 
by a minimal set of resonances and satisfy all constraints following from the heavy 
quark theory. This list can be continued. We are looking forward to new exciting 
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developments in the near future. 
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